in 
o 
o 

(N 



Open Quantum Dynamics: Complete Positivity and Entanglement 



OO 
(N 



Fabio Benatti 

Dipartimento di Fisica Teorica, Universita di Trieste, Strada Costiera 11 and 
Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 34014 Trieste, Italy; benatti@ts.infn.it 

Roberto Floreanini 

Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, Dipartimento di Fisica Teorica, 
Universita di Trieste, Strada Costiera 11, 34014 Trieste, Italy; florean@ts.infn.it 

fNj ' We review the standard treatment of open quantum systems in relation to quantum 

entanglement, analyzing, in particular, the behaviour of bipartite systems immersed in 
C D , a same environment. We first focus upon the notion of complete positivity, a physically 

l/^ i motivated algebraic constraint on the quantum dynamics, in relation to quantum en- 

C 2 3 ' tanglement, i.e. the existence of statistical correlations which can not be accounted for 

by classical probability. We then study the entanglement power of heat baths versus 
their decohering properties, a topic of increasing importance in the framework of the 
fast developing fields of quantum information, communication and computation. The 
presentation is self contained and, through several examples, it offers a detailed survey 
C ' of the physics and of the most relevant and used techniques relative to both quantum 

^ : open system dynamics and quantum entanglement. 
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1. Introduction 



5h 

' Standard quantum mechanics mainly deals with closed physical systems that can 

be considered isolated from any external environment, the latter being generically a 
larger system consisting of (infinitely) many degrees of freedom. The time-evolution 
of closed systems is described by one-parameter groups of unitary operators em- 
bodying the reversible character of the dynamics. On the contrary, when a system 
S interacts with an environment E in a, non-negligible way, it must be treated as 
an open quantum system, namely as a subsystem embedded within E, exchanging 
with it energy and entropy, and whose time-evolution is irreversible.^ 

In general, the time-evolution of S is inextricably linked to that of E. The 
compound system S* + £' is closed and develops reversibly in time; however, the 
global time-evolution rarely permits the extraction of a meaningful dynamics for the 



''The literature on the theory of open quantum systems and their phenomenological applications is 
vast; in the References, only papers that are strictly relevant to the present exposition are therefore 
included. We start by providing a list of general reviews and monographs on the topic, Refs.[l-16], 
where additional references can be found. 
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system S alone. This can be done if the couphng among subsystem and environment 
is sufficiently weak, in which case physically plausible approximations lead to reduced 
dynamics that involve only the S degrees of freedom and are generated by master 
equations. Such reduced dynamics provide effective descriptions of how E affects 
the time-evolution of S which, on time-scales that are specific of the given physical 
contexts, typically incorporates dissipative and noisy effects. 

In order to better appreciate the physical motivations underlying the concept 
of reduced dynamics, classical Brownian motion is a helpful guide. On the fast 
microscopic time-scale, the scattering of the environment particles off the meso- 
scopic Brownian particle S are described by a reversible dynamics. On the slower 
mesoscopic timc-scalc, the memory effects related to microscopic interactions are 
averaged out and the effective dynamics of S is diffusion, that is an irreversible 
Markov process. Physically speaking, on the slow time-scale the environment de- 
grees of freedom act as a sink where S dissipates energy, but also as a source of 
(white) noise with long-run stabilizing effects. We shall see that, under certain condi- 
tions, irreversible Markovian time-evolutions, leading to dissipation and noise, find 
a concrete description in quantum mechanics by means of the so-called quantum 
dynamical semigroups. 

Classical Brownian motion indicates that, when the typical time-scale of S is 
much larger than the time-scale governing the decay of time-correlations of the en- 
vironment, then E can be described as an effective source of damping and noise. 
In the framework of open quantum systems, this possibility is technically imple- 
mented either by letting the typical variation time of S, ts, go to infinity, while the 
environment correlation time te stays finite, or by letting te go to zero, while ts 
stays finite.'' As we shall see, these two regimes give rise to two different proccdiircs 
to arrive at a reduced dynamics described by Markovian semigroups: the so-called 
weak coupling and singular coupling limits. 

Since their first appearance, open quantum systems have been providing models 
of non-equilibrium quantum systems in diverse fields as chemical-physics, quantum 
optics and magnetic resonance. Recently, the rapid development of the theory of 
quantum information, communication and computation^ has revived the interest 
in open quantum systems in relations to their decohering properties, but also in 
their capacity of creating entanglement in multi-partite systems immersed in cer- 
tain environments. The typical open quantum systems in these contexts are n-level 
systems, like atoms, photons or neutrons embedded in optical cavities or heat baths 
consisting of bosonic or fermionic degrees of freedom. These will be the cases studied 
in this review. 

It is worth mentioning that, more in general, one can expect dissipative and 

noisy effects to affect also the behaviour of c!lc!ni(;ntary particle systems, like nciitral 
mesons, neutrinos and photons: the ensuing irreversibility is the emerging result of 
the dynamics of fundamental degrees of freedom at very short distances, typically 



''For general reviews on these topics, see Refs. [18-20]. 
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the Planck-scale.^^ These effects produce distinctive signatures in the behaviour 
of selected physical observables, allowing direct experimental testing (for specific 
treatments and further details, see Refs. [22-32]). 

Providing different concrete experimental contexts to study dissipative quantum 
dynamics is important not only from the point of view of the physics of open 
quantum systems, but also for the investigations of certain still debated aspects of 
these generalized time-evolutions, specifically rcigarding the property of positivity 
v.s. that of complete positivity, notions on which we now briefly elaborate. 

A primary consequence of a noisy environment E is that pure states of S, that 
is projections onto Hilbert space vectors, are usually mapped into mixtures of pro- 
jections, that is into generic density matrices. These are operators with discrete 
spectrum consisting of positive eigenvalues summing up to one; within the statisti- 
cal interpretation of quantum mechanics, they represent probabilities and are used 
to calculate the mean values of all physical observables pertaining to S. In order to 
be physically consistent, any reduced dynamics must thus correspond to maps on 
the space of states of S that preserve for all times the positivity of the spectrum. 
In other words, any physically consistent reduced dynamics must transform density 
matrices into density matrices at all (positive) times. 

The simplest way to arrive at dissipative time-evolutions is to construct them 
as solutions of suitable master equations of Liouville type: dtpt — ^[pt], where 
L operates linearly on the density matrices pt describing the state of S at time 
t. These are usually obtained by first tracing away the environment degrees of 
freedom, an operation that leads to a dynamical equation plagued in general by non- 
linearities, secular terms and memory eff'ects. Non-linearities and secular terms can 
be eliminated by requiring the initial states of S to be statistically decoupled from 
an initial reference state of E. On the other hand, memory effects are expected to 
be relevant on short time-scales and to become negligible if one looks at phenomena 
associated to longer time-scales when those effects have already died out. 

In technical terms, studying the time-evolution on a slow time-scale, thus ne- 
glecting memory effects, means operating a so-called Markovian approximation; this 
gives an efficient description of the reduced dynamics of S in terms of time-evolutions 
consisting of one-parameter semigroups of linear maps 7* = exp(tL), t > 0. These 
are maps that, acting on all possible initial density matrices p, provide their time- 
evolved partners pt = 7* [p] at time t. The dynamical maps 7t are formally obtained 
by exponentiating the generators L, so that 7t+s = ° 7s = 7s ° 7ti s,t>0. 

In most open systems, the memory effects are due to a transient regime with a 
very small time-span that makes Markov approximations legitimate; nevertheless, if 
not carefully performed, a likely outcome is a physically inconsistent time-evolution 
that does not preserve the positivity of the spectrum of all initial density matrices. 
As we shall explicitly see in the following, brute force Markov approximations may 
provide semigroups of linear maps 7f that at some time transform an initial state p 
into an operator pt with negative eigenvalues, which can not be used as probabilities. 
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In other words, without due care in the derivation, one is very hkely to obtain 7* 
that do not transform the space of states into itself. In order to avoid the appearance 
of negative probabihties, one has to perform Markov approximations that lead to 
semigroups of linear maps that preserve positivity: maps with such a property are 
called positive. 

Remarkably enough, both the mathematical property of linear maps of being 

positive and the derivation of positive dynamical maps through Markov approxi- 
mations are elusive and still not fully understood. However, positivity is not the 
end of the story and even a semigroup of positive linear maps would not be fully 
physical consistent. Indeed, suppose S is coupled to an inert system A, a so-called 
ancilla, with which it does not interact; then, the time-evolution of the compound 
system S + A is idA- The identity map id^ means that the coupling to the 
ancilla does not possess a dynamical character: it manifests itself purely in terms of 
statistical correlations, i.e. in the possibility of initial correlations between 5* and A. 

It turns out that the positivity of 7* does not guarantee the positivity of 7t(g)idA; 
namely, even if transforms any initial state of S into a state, 74 (8) id^ can map 
an initial state of S + A into an operator which exhibits negative eigenvalues and 
can not thus be interpreted as a physical state oi S + A. In order to be sure that, 
for whatever ancilla A, the evolution 74 (S) id^ consistently maps any state of -|- A 
into a state, 74 must be a so-called completely positive map for all t>0, a stronger 
property than positivity. 

We shall show that the ultimate physical meaning of complete positivity lies 
in its relation to the existence of entangled states, the foremost instance of them 
being a vector state with a singlet-like structure that can not be written as a tensor 
product of vector states. Indeed, the only states that could be mapped out of the 
state-space by jt id^ and thus made not anymore acceptable as physical states, 
are only the entangled states of S + A. 

Unitary timc-cvolutions are automatically completely positive, while this is not 
in general true of many Markovian approximations encountered in the literature 
on open quantum systems. Complete positivity or its absence depend on the way 
Markovian approximations are performed. Since the coupling of the open system 
S with a generic ancilla can not be excluded, the request of complete positivity 
is hardly dismissible when one deals with dynamical semigroups. Such a request 
has strong mathematical and physical consequences in that the structure of the 
generator L is fixed; from this a typical hierarchy among the characteristic relaxation 
times of 5* follows. The uncontrollability of the ancilla and the abstractness of the 
whole justification of why complete positivity should be physically compelling is 
scarcely appealing from a concrete point of view, all the more so since complete 
positivity strongly constrains the reduced dynamics. 

However, we shall see that the same physical inconsistencies typical of the ab- 
stract contexts, S + A, also occur when one considers not a generic inert ancilla A, 
but another system S also immersed in E, that is an open quantum system S + S. 
If the two parties weakly interact with E and do not interact between themselves. 
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the dissipative time-evolution of the compound system 5 + 5 is of the form 7* (g) 7^ 
instead of the abstract jt •8 id^ and the setting becomes more physical. 

Consider, for instance, two optically active non-interacting atoms (or molecules) 
initially prepared in an entangled polarization state and evolving in contact with a 
same heat bath, weakly coupled to them. Then, if the dynamical maps 7* describing 
the dissipative evolution of each single atom is not completely positive, there are 
physically admissible initial entangled states of the two atoms that, after some time, 
lose their interpretation as physical states because of the emergence of negative 
probabilities generated by ^t®lt- 

The rising of quantum entanglement as a physical resource enabling the per- 
formance of quantum information and computation tasks otherwise impossible as 
teleportation, superdense coding, quantum cryptography and quantum computa- 
tion, has spurred the investigation of how entanglement can be generated, detected 
and manipulated. Contrary to standard expectations, the presence of an environ- 
ment, and thus of noise, need not necessarily spoil the entanglement properties of 
states of systems immersed in it, but in some cases can even have entangling effects. 
Using the previous examples, the two optically active atoms could be prepared in 
an initial separable state without either classical or quantum correlations and put 
into a heat bath. The covariance structure of the bath correlation functions may be 
such that not only the two atoms become entangled, but also such that a certain 
amount of entanglement survives over longer and longer times. The reason is not 
hard to see: even if not directly coupled, two systems immersed in a same heat bath 
can interact through the bath itself; it does then depend on how strong this indirect 
interaction is with respect to the decoherence whether entanglement can be created 
and maintained. 

The plan of this review is as follows. In Section 2, we shall set the mathematical 
framework appropriate to the description of open quantum systems and to the 
derivation of their reduced dynamics, with particular emphasis on the notions of 
positive and completely positive maps on the space of states. 

In Section 3, we shall review standard derivations of dissipative semigroups, 
discussing the role of the Markov approximations in achieving completely positive 
dynamical maps or not and illustrating the physical inconsistencies arising from 
lack of complete positivity. 

In Section 4, we shall address the entangling effects that the presence of an 
environment may have on non-directly interacting bipartite systems. All topics will 
be illustrated by means of concrete examples. 

2. Open Quantum Systems: Mathematical Setting 

In the first part of this Section we review some basic kinematics and dynamics of 
open quantum systems, fixing, in passing, the necessary notation. In the second part 
we shall discuss positivity and complete positivity of linear maps with particular 
emphasis on their relations with quantum entanglement. 
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We shall consider finite (n- level) open quantum systems S; they are described by 
means of n-dimensional Hilbert spaces C", where C is the set of complex numbers, 
and by the algebras M„(C) of n x n complex matrices X.'^ 

The hcrmitian matrices X = X'^ correspond to the system observablcs and their 
mean values depend on the physical states of S. The latter can be divided into two 
classes, the pure states described by projectors = IV') (V"! onto normalized vectors 
Irp) € C" and the statisfAcal mixtures described by density matrices, i.e. by linear 
convex combinations of (not necessarily orthogonal) projectors 

p = ^A,P^. , A,>0, ^A, = l. (1) 

3 3 

In the following we shall refer to them simply as states, specifying whether they are 
pure or mixtures when necessary; further, 

TV(p) = $^A, = l, (2) 

3 

where Tr represents the trace-operation Tr{X) = 'Y^^{i\ X \i), with i = 1,2, ■ ■ ■ ,n 
any orthonormal basis in C". 

Projectors and density matrices are thus normalized hermitian elements of 
M„(C); further, they are positive semi-definite matrices. 

Definition 2.1 A matrix X^ = X G M„(C) is positive semi-definite if 

(V'l X IV-) > , V V e C" . (3) 

Positive semi-definite matrices will be called positive for sake of simplicity and 
denoted by X > 0; their spectrum necessarily consists of positive (non-negative) 
eigenvalues. 

Quantum states are thus positive, normalized operators; the eigenvalues of the 
pure ones are 1 (non-degenerate) and (n — 1 times degenerate), while those 
of density matrices are generic < rj < 1, j = 1,2, ... ,n, such that J2j ''i = 1- 
It follows that pure states and mixtures are distinguished by the fact that = 
p if and only if p = for some |V') G C". Also, the eigenvalues of a density 
matrix p correspond to the weights Xj in (1) if and only if (1) is the spectral 
representation of p, the Itjjj) are its eigenstates and the P.^. the corresponding 
orthogonal eigenprojectors, P^-P^,^ = SjkP^^- 

Given a state p, the mean value {X)p of any observable X = e M„(C) is 
calculated as follows: 

(X),^Tr(pX)=^A,(V^,|X|V,) . (4) 

3 

The previous considerations constitute the bulk of the statistical interpretation 
of quantum mechanics: 



''General monographs on quantum physics relevant for our exposition are Refs. [33,36,37]. 
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Given the spectral decomposition of a state, 

n n 

p = "^ri \ri){ri\ , = 1 , {ri\rj) = % , (5) 

i=l i=l 

the eigenvalues constitute a probability distribution which completely defines the 
statistical properties of the system. 

Accordingly, one associates to quantum states the von Neumann entropy, 

n 

S{p) = -Tr(plogp) = - ^ r, logn , (6) 

which measures the amount of uncertainty about the actual state of S. It turns out 
that S{p) = if and only if = p, otherwise < S{p) < logn. 

Remark 2.1 Density matrices form a convex subset S{S) c M„(C) which we 
shall refer to as the state-space of S. Namely, combining different mixtures aj e 
15(5) with weights \j > 0, = 1, into the convex combination ^j<^j, the 

latter also belongs to S{S). Pure states are extremal elements of S{S), that is they 
can not be convexly decomposed, while with them, by linear convex combinations, 
one generates the whole of the state-space. □ 



2.1. Reversible and Irreversible Dynamics 

The state-space and the algebra of observables fix the kinematics of S. Its dynamics 
as a closed system is determined by a Hamiltonian operator H e M„(C) through 
the Schrodinger equation (we shall set H= 1): 

dt\i,t) = -iH\i,t) . (7) 

By direct inspection, first on projectors and then on mixtures, this gives rise to 
the so-called Liouville-von Neumann equation on the state-space S{S): 

dtpt = -I [H, pt] , (8) 

whose solution, with initial condition pt=o = p, is 

pt = UtpU.t, Ut = e-'"'. (9) 

Denoting by p U( [p] = pt the dynamical map (9) and by 

p^-LH{p] = -i[H,p] , (10) 

the linear action of the generator on the left hand side of (8), the Schrodinger unitary 
dynamics amounts to exponentiation of L: 

j-k 

Pt = Ut[X] = e*"^^ [p] = V - Lh o • • • o LhM , (11) 

ki •• V ' 
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where o means compositions of maps. Therefore, the dynamical maps Uj form a 
one-parameter group of linear maps on S{S): Ut o Ug = Vt+s for all t,s gR. This 
fact mathematically describes the reversible character of the unitary Schrodinger 
dynamics; in particular, the dynamical maps Ut can be inverted, preserve the spec- 
trum of all states p, leave the von Neumann entropy invariant and transform pure 
states into pure states: 

p' = p=^{Ut[p])' = Vt[p]- (12) 



Remark 2.2 One can pass from the Schrodinger to the Heisenberg time-evolution 
through the definition of mean values (4) and the so-called duality relation 



Tr 



= Tr 



pVliX] 



(13) 



which holds for all p G S{S), X e M„(C) and defines the dual linear map U^. The 
latter acts on M„(C) as 



X ^ U* [X] = U-t XUt = e-*'^« [X] 



(14) 
□ 



The unitary dynamics is not the only dynamical transformation affecting quan- 
tum states. According to the postulates of quantum mechanics, if the state of S is 
\ip){'tjj\, upon measuring the (spectralized) observable X = Xi\xi){xi\, then 

• the eigenvalues obtained with probabilities |(a;,|'^)|^, 

• if the measure gives Xi, then the post-measurement state of S is |a;8)(a;j|. 

It follows that, by repeating the measurement many times on copies of S equally 
prepared in the pure state |V')(V-'I ^.nd collecting all the resulting post-measurement 
states, the outcome is a physical mixture described by the density matrix 

n n 
i=l i=l 

Setting Pi = \xi){xi\ and extending (15) linearly to any density matrix p G S{S), 
one gets the following linear map on S{S): 

n 

p^V[p]=Y,PipPi ■ (16) 

i=l 

This map is a mathematical description of the so-called wave-packet reduction. 

Contrary to the unitary dynamics Vt, P transforms pure states into mixtures 
thus increasing their von Neumann entropy; the process it describes is sometimes 
identified with a randomizing quantum noise, with decohering properties. By deco- 
herence it is meant the result of any dynamical transformation that suppresses the 
phase-interferences present in a linear superposition of vector states. 
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Remark 2.3 The wave-packet reduction mechanism P effectively describes what 
happens to S when it is not closed but in interaction with an external specific 
environment, in this case an apparatus measuring the observable X. Specifically, 
the wave packet reduction is the final effect on the open system S of its interaction 
with the macroscopic environment. □ 

In the following, we shall be concerned with open quantum systems S immersed 
in an environment E. In principle, the environment should consist of infinitely many 
degrees of freedom and thus be addressed by means of the more proper algebraic 
approach to quantum statistical mechanics;'^^'^^ however, for sake of clarity, we shall 
describe it by means of an infinite dimensional Hilbert space Sj and represent its 
states by density matrices p^.^ 

Subsystem S and environment E make a closed compound system S + E, its 
Hilbert space being the tensor product C" (Xi ^5; the total system evolves revcrsibly 
according to a group of dynamical maps Uf^^ = ex.p{tL,s+E) that act on the state- 
space S{S + E). The group is generated by formally exponentiating the commutator 
^s+e[ps+e] = —i[Hs+E , Ps+e] with respect to a total Hamiltonian 

Hs+e = Hs<»1e + 1s<»He + XH' , (17) 

where A is an adimensional coupling constant, H' describes the S — E interaction, 
while Hs and He are Hamiltonian operators pertaining to S*, respectively E and 
Is.B are identity operators. It follows that the generator decomposes as a sum of 
commutators, 

L = Ls+B =Ls Lb -h AL' , (18) 

where the subscripts S, E identify which degrees of freedom are involved. 

Given a state ps+E G S{S + E) of the compoimd closed system, the statistical 
properties of the embedded subsystem S are described by a state ps G S{S) which 
is obtained by the partial trace over the degrees of freedom of E: 

S{S + E)3 ps+E ^Ps = Tte(ps+e) = J2 (^f I Ps+E IV'f ) , (19) 

3 

where {IV'j')} is any orthonormal basis in 9). The right hand side of (19) belongs 
to M„ (C) and can be easily checked to be positive and normalized according to (3) 
and (2), so that ps € S{S). 

Analogously, given any ps at time t = 0, the state of S at any time t is 

ps{t) = TrE{uf+''[ps+E]) ■ (20) 
This gives rise to a family of maps Gt, 

PS ^ Ps{t) = Gt[ps] , (21) 



At a certain stage, one nevertheless has to perform the so-called thermodynamic limit, in which 
volume and number of degrees of freedom are let to infinity, while keeping the density finite. In prac- 
tical terms, this amounts to substituting integrals for discrete summations, e.g. see Example 3.3. 
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which in general depend on ps and can not be extended to the whole of the state- 
space. If we ask that preserve the convex structure of S{S), that is 

Gt[EVs] =E^^-^*k^] ' ^^^^ 

j 3 

then the initial state of the compound system must factorize: ps+E = Ps ® Pe, 
where ps G 5(5') and pE is a fixed state of the environment.''^ 

The factorized form of the initial state means in particular that open system 
and environment are initially completely uncorrelated; though not true in general, 
in many interesting physical contexts such a condition is fully consistent and gives 
rise to a family of dynamical maps Gt depending on the environment reference state 
Pe, but otherwise acting linearly on the state-space of S.^ 

The family of maps Gt, t>0, describes a forward-in-time irreversible dynamics, 
for the partial trace breaks time-reversal; what the family lacks is a semigroup 
composition law, since in general G* o Gg ^ Gt+s for t,s > 0. An equality would 
express the absence of cumulative memory effects, as such it is expected to be 
a good approximation of the time-evolution of open quantum systems only when 
their interaction with the environment is sufficiently weak or the environment time- 
correlations decay rapidly with respect to time- variation of the system. 

The technical procedures to eliminate memory effects and recover semigroups of 
dynamical maps as reduced time-evolutions are known as Markov approximatAons. 
These will be discussed in detail in the next Section; in the remaining part of the 
present Section we shall study the structure of Gt as a linear map, focusing on its 
physical properties. 



2.2. Positivity and Complete Positivity 

As motivated before, we take an initial state of the compound system S + E in 
factorized form ps+E = Ps ® Pe- By performing the partial trace with respect to 
the orthonormal basis of eigenvectors \rf) of pE with corresponding eigenvalues rf , 
one obtains 

ps{t) = IV^(uf+^[ps ® pe]) = Y: r^{rf\Ur^\r^ ) Ps (rf |C/f,+^|rf ) 

3,k 

= J2Vc.{t)psVUt) , (23) 

a 

where a is a double summation index and 

Vo.it) ^ ^ {rf\Uf+^\ri} . (24) 



''As already mentioned, in presence of initial correlations between subsystem and environment, a 
reduced dynamics that acts linearly on all S{S) can not in general be defined: its explicit form 
crucially depends also on the initial reduced state TteIps+e]', for further investigations on this 
issue, see Refs.[40, 41, 7], and references therein. 
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The matrix elements {'i'f\Uf^^\rf) are elements of M„(C); further, from normal- 
ization (Tr(ps(i)) = 1), it follows that 

Y.Vi{t)Vo.{t) = ln, (25) 

a 

where 1„ e M„(C) denotes the identity matrix. 

The duality relation (13) associates to the linear map 

PS ^ [ps] = ^ Va {t) psVlit), (26) 

a 

on the state-space S{S)^ the dual action 

X ^ G* [X] =J2vi (t) X Va it) , (27) 

a 

on the algebra of observables Af^((C). The latter fulfils (C-^[ly7,] = ly^, a property 
called unitality. 

The structure of both and G^ is similar to the wave-packet reduction (16), 

with the difference that the Va{t) need not be orthogonal one-dimensional projectors 
as the Pi. By means of (3), it is easy to see that these linear maps preserve the 
positivity of operators: they belong to the class of positive maps.^^ 

Definition 2.2 A linear map A : M„(C) M„(C) is termed positive if it sends 
positive matrices into positive matrices, namely if 

0<X ^ A[X] > . (28) 



Remark 2.4 While positivity of hermitian matrices means positivity of their 
eigenvalues, positivity of linear maps means their property of transforming a ma- 
trix with positive eigenvalues into another such matrix. Positive linear maps are 
sometimes more properly referred to as positivity-preserving maps. □ 

It turns out that linear maps as Gf, G^ and P belong to a special subclass of 
positive maps on M„(C): the so-called completely positive maps.^^'^^ 

Complete positivity is a stronger property than positivity. It concerns the possi- 
bility that the system S be statistically coupled to a so-called ancilla, A, that is to 
a generic remote and inert, finite-dimensional system. If A is m-dimensional, any 
linear map A on Af„ (C) lifts to a map Aj^ id^ on the algebra M„(C)(8)Mm(C) of the 
compound system S + A, whore id^ denotes the identical action on M„i(C) leaving 
all operators unaffected. In the following the identical action will appear either as 
id^ or as id^, explicitly indicating the dimensionality of the ancillary system. 

The physical interpretation of such a coupling is as follows: the system S of 
interest may have interacted in the past with A and become correlated (entangled) 
with it; afterwards and prior to the start of the evolution of S, S and A has ceased 
to interact so that what is left are the statistical correlations between them incorpo- 
rated in the state ps+A- When the evolution (or more in general any state change) 



12 Benatti, Floreanini 



of S sets in, it has to be considered as a transformation of the compound system 
S + A, which does not affect the ancilla and is thus mathematically described by 
the lifting A (g) id^ of the map A. 

Definition 2.3 A linear map A : M„(C) i-^ M„(C) is completely positive, if and 
only if id^ is positive on M„(C) (g) Mm{C) for all m > 1, that is for all possible 
statistical couplings with finite- dimensional ancillas. 

When m = l, 5 + ^is5 itself as the Hilbert space C" ® C is isomorphic 

to C". When m > 2, all elements of M„(C) (E) Mm(C) can be written as m x m 
M„(C)-valucd matrices Xij £ M„(C); therefore, the definition concerns 

the linear actions A (g) idA[[Arjj]] = [A[Xjj]] and selects among the positive linear 
maps A those such that*^ 





/Xn ■ 
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In the following Theorem, wc recall two important results characterizing com- 
pletely positive maps (sec Refs.[42, 34, 33] for their proofs). The first one states 
that, in order to ascertain the complete positivity of A : M„(C) ^ M„(C), it is not 
necessary to check all dimensions m, while the second one establishes that complete 
positive A are exactly those with the same structure as Ut, P, and G^. 

Theorem 2.1 

• A linear map A : M„(C) i— > M„(C) is completely positive, if and only if 

A (g) id„ is positive on Mn (C) (g M„ (C) . 

• A linear map A : M„(C) ^ M„(C) is completely positive, if and only if it 
is expressible in Kraus-Stinespring form 

Am=E^«^"^«' (29) 

a 

where Va € M„(C) are such that J2a ^a^a converges (to 1„ if A is unital). 

Suppose one has a linear map on M„(C) given in the form 

A[X] = J2Ca0WiXW0 , (30) 

with Wq, g M„(C) and Cafj making a hermitian matrix of coefficients such that 
J2ai3 ^af3 converges. By diagonaUzing [Cafj] = U^dmg{di.d2. . . .)U, one re- 

covers (29) if and only if [Ca/s] is positive definite and the eigenvalues di positive; 
indeed, Cq,/? = J2j djU^^Uj/j and this serves to redefine Vj = ^/d]Uj^yVa■ 
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Remark 2.5 If not all dj are positive, then the map can be separated into two 
sums, one over the positive and the other over the negative eigenvalues; by extracting 
an overall minus sign. A can then be written as the difference of two completely 
positive maps A1.2: A = Ai — A2. The property of complete positivity fixes the 
structure of the map by excluding the presence of a A2, while that of positivity is 
still far from being understood: no general prescriptions on Ai^2 are known that 
ensure the positivity of A. □ 

Completely positive maps are a proper subset of all positive maps, for not all positive 
maps A lift to positive A (g) id^, the foremost example being the transposition. 



Examples 2.1 

(1) In M2(C), let T2 : 



a b 
c d 



b d 



; the transposition map T2 does not 



alter the spectrum and therefore is a positive linear map. Let |0) 

'0' 



|1) 



be an orthonormal basis in and consider the vector 



|4.(^2)) = -L(|0)(^|0) + |l)(g|l)) eC2(^c2 . 



p(2) 



1)(1|®|1)(1| + |0)(1|®|0)(1| + |1)(0|®|1)(0| 



The corresponding projection = |*^^)(*^^| e M2(C) O M2(C), 

1 
2 



/I 1\ 





\l 1/ 

has eigenvalues 0, 1 and gets transformed by T2 ® id2 into 



(31) 



T2®id2[Pf] = i(|0)(0| 



+ |1)(1|«|1>(1| 

+ |1)(0|®|0)(1| + |0)(1|®|1)(0|) 



/I 0\ 

10 

10 

Vo 1/ 



which has eigenvalues ±1/2 and is no longer a positive matrix. Therefore, T2 is 
positive, but not completely positive. 
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(2) The vector state |^'^ ) generalizes to the totally symmetric state 

|j)®|j)gC"0C" , (32) 

where \j), j — 1,2, ... ,n, is a. reference orthonormal basis in C". Then, with 
p|") = |*^"^)(*^"^|, and T„ the transposition on M„(C), 

and one checks that the sum in the right hand side of the second equality is the 
flip operator V onC"(g)C", such that V (S) {(j)}) = {(j)} ® , W \^} , {(j)} e C". 
It follows that V has eigenvalue —1 on any anti-symmetric state of C" (g) C" 
and thus the transposition is never completely positive in any dimension. □ 



Remark 2.6 Transposition is associated with time-reversal;^' as its lifting T„ 

(n) 

id„ sends the state P]_ out of the space of states, it can not correspond to a physical 
operation. This means that time-reversal has to be applied globally and not locally, 
that is not on subsystems alone. □ 



Examples 2.2 Let 5 be a quantum system described by a 2-dimensional Hilbert 
space C^; the algebra of observables M2(C) is linearly spanned by the identity 
matrix and the three Pauli matrices 















I) 









-i 



0-3 



1 
-1 



namely, X = X]^=o '^fJ- '^ith ^/i — 5Tr(X(7^). The linear map 
Sa : M2(C) M2(C) , X §o,[X] = a^X aa , 

are written in the Kraus-Stinespring form (29), and are thus completely positive 

(although not unital). 



(1) Consider the following two maps, 



3 3 



(33) 



where £„ 



Q = 

1 when aj^ 2, whereas £2 



a=0 



-1. Using the algebraic relations 



aa cr/3 ctq = r?a/3 , where = 1 



f3\a 





1 


2 


3 





1 


1 


1 


1 


1 


1 


1 


-1 


-1 


2 


1 


-1 


1 


-1 


3 


1 


-1 


-1 


1 
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it can be checked that the first map amounts to X ^ Tr{X) ctq and is com- 
pletely positive since it is in Kraus-Stinespring form, while the second one is 
the transposition which indeed is such that T2[(Ta] = £a<^a, therefore positive, 
but not completely positive. 
(2) Given A : M2(C) ^ M2(C) in the form (30), by expanding the along the 
Pauli matrices aa it follows that the map can be rewritten as 

3 

X ^ A[X] = <T^Xa0, (34) 

a, 13=0 

and is completely positive if and only if the matrix [Ca^] is positive definite. □ 

The different structure of positive and completely positive maps on M„(C) is better 
exposed by means of the state Pj."'' of Example 2.1.2; one has:''^ 

Theorem 2.2 

• A : M„(C) I— > M„(C) is positive if and only if 

(V'0?!'|A(g)id„[P|"^]|VO?!') >0 , VV',«!)eC". (35) 

• A : M„(C) M„(C) is completely positive if and only if 

(*| A id„[P|"^] I*) > , V* e C" O C" . (36) 

The request in (36) is that the matrix A0 id„[_p|"^] G M„2(C) be positive (see 
Definition 2.1), while the request in (35) involves mean values with respect to prod- 
uct vector states, only. These latter do not exhaust the Hilbert space, as there are 
non-product states as 

2.3. Complete Positivity and Quantum Entanglement 

The projector pj"-' is the typical instance of an entangled pure state, the general 
notion of entanglement as opposed to that of separability being as follows. 

Definition 2.4 Let S{Si) and S{S2) be the state-spaces of two n-level quantum 
systems S'1,2; within the state-space S {Si +S2) one distinguishes the (convex) subset 
of separable states of the form 

= J2 ^ p^i ' > , ^ = 1 , p\:] e 5(^1,2) . (37) 

All those states which can not be written as linear convex combinations of product 
states as in (37) are called entangled. 
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Remarks 2.7 

(1) The first two terms in the first line of (31) represent separable states: it is the 
presence of the remaining two interference terms that makes Pj.^^ entangled. 

(2) It is often stated that while separable states only carry the classical correlations 

embodied in the weights Xij in (37) , entangled states carry instead non-classical 
correlations, that is they arc characterized by statistical properties which are 
not describable by classical probability theory. This is best seen by considering 
the projector p|"^ = '|, its von Neumann entropy S(P|i"'') = 0. The 

(n) 

partial trace of P| over either the first or the second system in S* + 5* gets the 
totally mixed state t„ = In/n. This state has maximal entropy 5'(t„) = logn, 
despite S being a subsytem of S+S, which has zero entropy; instead, in classical 
probability theory, knowledge of the whole means having full information on all 
its parts. □ 

At first glance one may think that the positivity of linear maps A on S{S) may 
be sufficient to guarantee a consistent description of all possible physical transfor- 
mations of the states of S: after all, if a map A preserves the trace and the positivity 
of the spectrum of all density matrices, then it maps the state-space S{S) into itself 
as physical transformations surely do. In fact, a minimal consistency request is 

Physically consistent transfornnations on S{S) must preserve the interpretation of 
the eigenvalues of all density matrices as probabilities and thus respect the positivity 
of their spectra. 

One can now appreciate the important physical consequence of Theorem 2.2; 
indeed, from (35) and (36) it turns out that if A : M„(C) M„(C) is positive, but 
not completely positive, then A (g) id„[Pj"^] is not a positive matrix. 

Given an n-dimensional ancilla A, statistically coupled to S, operating the trans- 
formation A on S{S) amounts to act with A <S) idA on the state-space S{S + A) of 
the compound system S + A. Now, among the states oi S + A, the totally symmetric 
projector pj"^ is a physically plausible initial state, which is not transformed into 
a state by A (g) id^- For this to be true, A has to be completely positive. Positiv- 
ity alone suffices for physical consistency on S, but not on S + A; in other words, 
positive maps can not properly describe physical transformations, only completely 
positive maps may do. Hence, 

Complete positivity of maps is necessary to ensure their physical consistency against 
possible entanglement with ancillas 

Remark 2.8 Admittedly, though logically stringent, the explanation of why con- 
sistent state-transformations must be completely positive is scarcely appealing from 
a physical point of view. Indeed, the ancilla system A is remote from S and inert; a 
more concrete scenario will be offered in Section 3.2 when we deal with dissipative 
time-evolutions. On the contrary, couplings to ancillas are natural tools in quantum 
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information, although, in this case, one deals with quantum channels and not with 
continuous time-evolutions. □ 

If positive maps do not consistently describe physical transformations, they play 
however a fundamental role in the detection of entangled states. Due to the re- 
cent fast advances in quantum information theory, entanglement has turned from 
a quantum riddle with cpistcmological overtones into a practical physical resource 
that must be detected, quantified and manipulated. For instance, bipartite entan- 
gled pure states as P]_ ' are essential ingredients in quantum cryptography, quantum 
teleportation and quantum computation.^^ 

2.3.1. Entanglem,ent Detection 

From Definition 2.4, entangled states can generically be mixed. While, it is easy to 
spot entangled pure states, it is much harder in the case of density matrices. 

Example 2.3 As a first simple indication of these difficulties, consider the totally 
depolarized state of two 2-dimensional systems T2 = I4/4. The identity matrix can 
be written in terms of the projectors onto the orthonormal basis of the so-called 
Bell states, 



where Pf^ is as in Example 2.1.1 while P^^ = Q^^^ and 

I*?)) = 4(100) -111)) , |4^)) = 4(|oi)±|io)) . 



As a combination of entangled projectors, T2 may be deemed entangled. However, 
mixtures of entangled states have less entanglement than their constituents; indeed. 



and thus, as a combination of separable states, T2 is indeed separable. □ 

A substantial help in detecting entanglement comes from the transposition as 
defined in Examples 2.1, where one looked at the effects of the action of the trans- 
position on only one of the factors of a bipartite system S + S, the so-called partial 
transposition T„ (gi id„. Also, based on the fact that separable states make a closed 
convex subset of the state-space which can be geometrically separated from any 
given entangled state by a suitable (hyper) plane, one deduces the so-called Peres- 
Horodecki criterion for separability.^'^ 

Theorem 2.3 Let S he an n-dimensional system. Then, 

• a state p G S{S + S) is entangled if it does not remain positive under partial 
transposition T„ id„.^'* 
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More in general: 

• a state p e S{S + S) is entangled only i/ 3 A : M„(C) M„(C) positive such 
that A (g) id,i [p] is not positive.^^ 

The sufficient condition easily follows from the fact that the action of any positive 
map A, lifted on separable states psep as in (37), 

A O id„ \psep] = ^ Xij A[pI] O p'j , (38) 

keeps their positivity, for A keeps the positivity of all p^ G »5(S'i) by assumption. 
Therefore, if a positive map like the transposition, is such that once lifted it does 
not preserve the positivity of p e S{S + S) then the state must be entangled: such 
a positive map is called an entanglement detector for p. 

Of course, entanglement detectors can not be completely positive maps since, 
according to Definition 2.3, they are exactly those positive maps on M„(C) such 
that, once lifted to act on M„(C) (g) M„(C), they remain positive. 

In line of principle, the transposition is a good entanglement detector for some 
states and not for others; one may thus suspect that, in general, all positive maps 
have to be checked in order to detect the entanglement of a generic p. Fortunately, 
in the case of a bipartite system consisting of two 2-lcvcl systems, or more in general 
of one 2-level and one 3-level system, the property of not remaining positive under 
partial transposition is also a necessary condition for their states to be entangled. 
In other words,^^ 

Theorem 2.4 Consider a bipartite system Si + S2, with Si a 2-level system and 
S2 either a 2-level or a 3-level system; a state p G S{Si + S2) is entangled if and 
only iff 2 (g id2[p], respectively T2 (g idafp], is not positive. 

ExEimple 2.4 Consider the x matrix: 

+ /jy eM„(C)«)M„(C) , a,/?eR, 
where V is the flip operator introduced in Example 2.1.2, 

y = nT„(g)id„[P|"^] , (39) 

and the parameter F is defined by 

F = Tt[pfV] ^ an + 0ni^ . (40) 

The second equality follows from the fact that = 1„ (g 1„ and that, by choosing 
a basis \i,j) — \i) |j) e C" (g) C", then V\i,j) = We want pp to be a 

density matrix. Since V has eigenvalues ±.1, pp has eigenvalues a ± /3, positivity 
and normalization are guaranteed by 

-a</3<a, Tr(pir) =Q;n^ + /3n = 1 . (41) 
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Conditions (40) and (41) give 



n(n^ — 1) n[n^ — 1) 



I L ± I C J. ± ^ . 

PF = TT In «) In + -r^ 7T 1^ , 42 



The states of the one-parameter family 

n — F n F — 1 

n(n^ — 1) n{n'' — 1) 

are called Werner states.'^'^ They are all and the only states that commute with all 
unitary transformation on C" C" of the form U ®U with U any unitary operator 
on C". By performing the partial transposition on pF, using (39) one gets 

^ ,r, n — F nF — 1 r^) 

T„ ® idnlpp] = r 2 1^ In «> In + , 2 Pi ' , 
n{n^ — 1) (n^ — 1) ^ 

n — 

with an (n^ — l)-degenerate eigenvalue , „ which is never negative, and a 

n(n^ — 1) 

nondegencrate eigenvalue 

n- F nF -1 F 

+ 



n{rfi — 1) {v? — 1) n ' 

which is negative when F < 0. Therefore, in agreement with the preceding discus- 
sion, when F < 0, the Werner states are necessarily entangled, otherwise the partial 
transposition could not spoil the positivity of their spectrum. 

It turns out that F < is not only sufficient, but also necessary for to be 
entangled:^^ in other words, pi? is separable if and only if F > 0. Since F > means 
remaining positive under partial transposition, the latter is actually an exhaustive 
entanglement-detector relative to the family of Werner states also in dimension 
larger than two. □ 



In the case of bipartite 2-dimensional systems, there exists a quantitative method 
to determine whether a state p G S{S + S) is entangled or not, based on the notion 
of concurrence^^''^'^-^^ This is defined as follows: with the tensor products of Pauli 
matrices (72 ® cr2, one constructs 

ii = p (72 (72 p* (72 (72 , (43) 

where p* denotes p with complex-conjugated entries. The matrix R turns out to 
have positive eigenvalues; let Ri, i = 1,2,3,4, be the positive square roots of its 
eigenvalues in decreasing order. 



Theorem 2.5 A state p € S{S -\- S), with S two-dimensional, is entangled if and 

only if its concurrence 

C{p) = max|i?i - R2 - Rz - R4, o| , (44) 

is strictly positive. 
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Example 2.5 We consider the family of Werner states studied in the previous 
Example and set n = 2. Using the Bell states introduced in Example 2.3, the flip 
operator explicitly reads 



(2) 



)(2) _ 



whence 



PF = 







V 



/I 
















1 








1 








Vo 








1/ 







\ 




1 





















2 - P 2P ■ 

2P - 1 2 - 

1 + P/ 

Further, one cheks that a2 ® (J2 P f (^2 ® (^2 = Pf', then (43) gives R = p%; therefore, 
the eigenvalues Ri are (1 + P)/6 (three times degenerate) and (1 — P)/2, which 
are non-negative, since — 1 < P < 1; thus, the difference Ri — R2 — R3 — Ri 
in (44) can assume only two expressions: — P and (P — 2)/3. As a consequence, the 
concurrence is C{pf) = — P > 0, if and only if P < 0, which, as already seen in the 
previous Example, is a necessary and sufficient condition for pp to be entangled. □ 



3. Reduced DynEimics and Meirkov Approximations 

According to the previous considerations, the completely positive maps Gt in (26) 
constitute a physically consistent description of the irreversible dynamics of a system 
S in interaction with its environment E. The only necessary assumption is that the 
initial state of 5 + P do not carry either classical or quantum correlations and be 
of the form ps® Pe- 

However, the dissipative and noisy effects due to E are hidden within the opera- 
tors Va{t) and quite difficult to be read off; in particular, the one-parameter family 
of maps Gt contains memory effects. On the other hand, if the interaction between 
S and E is sufficiently weak, one expects that, on a typical time-scale, the dynamics 
of S might be disentangled from that of the total system and efficiently described 
by a one-parameter semigroup of maps 74, t > 0, satisfying the forward in time 
composition law 7( o 7, = 7t+s, t, s > 0. 

The memory effects present in are best revealed by writing the formal integro- 
differential evolution equation of which the maps Gt are solutions. Its standard 
derivation is via the so-called projection technique' which we shall review in some 
detail by presenting and comparing several Markov approximations, leading to mas- 
ter equations of the form 

dtpsit) = (Lff + mpsit)] , (45) 



'Originally introduced in Refs.[51, 52], it is widely used in non-equilibrium statistical mechanics. 
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where Ljj acts on the state-space S{S) as {—i times) the commutator with an 
effective Hamiltonian H = S M„(C) as in (10), whereas D is a linear operator on 
S{S), not in the form of a commutator, that effectively accounts for the dissipative 
and noisy effects due to E. 

The solutions 7t, t > 0, to (45) will describe the reduced dynamics of S" as a 
semigroup of linear maps on S{S), obtained by exponentiating the generator: 



7t = exp 



t {TLh + D) 



(46) 



The arguments developed in the previous sections about positivity and complete 
positivity lead to the conclusion that the latter property is necessary to ensure 
physical consistency against couplings to ancillas. Therefore, one can state in full 

generality that 

any physically consistent Markovian approximation must yield semigroups jt con- 
sisting of completely positive maps. 

In the following, we shall not abide by this request and expose the inconsisten- 
cies thereby arising; before this, however, we discuss the general form that time- 
continuous semigroups of completely positive maps must have. 

3.1. Abstract Form of Generators of Dynamical Semigroups 

Semigroups of the form (46) fulfill limt^o It = id; more precisely 

lim||7([p]-p||i = 0,VpG5(5) , 

where = TrVXtX, X G M„(C). 

Time-continuity of a semigroup of maps 74 guarantees the existence of a gener- 
ator and an exponential structure as in (46).'^^ We now proceed by imposing one 
by one those requests that are deemed necessary to physical consistency.^'^'^ 

The first constraint is that the hcrmiticity of density matrices be preserved; we 
shall also ask for trace preservation. Physically speaking, this means that the overall 
probability is constant; in other words, we shall not be concerned with phenomena 
like particle decays that are characterized by loss of probability.^ 

Interestingly, hermiticity and probability preservation sufHce to partially fix the 
form of the generator. 



Remark 3.1 Instead of directly referring to transformations on the state-space, 
the following theorem, as present in the literature, deals with semigroups 7^ acting 



SThese can nevertheless be included in the formalism: for instance, see Refs.[l, 25] and 
Remark 3.3.3 below. 
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on the algebra of observables A'/„(C) (* denotes dual transformations on observ- 
ables, namely in the Heisenberg picture). Hermiticity and trace preservation cor- 
respond to unitality, 7j*[l] = 1, and (it [X])^ = 7^ [Xt] for aU X e M„(C). Time- 
continuity is given with respect to the norm- topology: limt^o ||7t [-^] — -^|| = 0, for 
all X e M„(C), where 

= sup{^(vixtxiv.) , \m = 1 } • 

Upon passing from the Heisenberg to the Schrodingcr picture, one gets semigroups 
acting on states by means of the duality relation (13). □ 

Theorem 3.1 Let 7^* : M„(C) 1— > M„(C), t > 0, form a time- continuous semi- 
group of unital, hermiticity-preserving linear maps. Then, the semigroup has the 
form 7t = exp(f(L|^ + P*)) with generator consisting of 

h%[X] = i[H , X] , (47) 

— 1 ^ 1 \ 

= J2 Cij ( Fi X f] - -{PiF} ,X} \ , (48) 
i,j=i ^ ^ 

where the matrix of coefficient Cij (Kossakowski matrix) is hermitian, H = H\ 
and the Fj are such that F„2 — In/y/n and Tr{Fj Fk) = Sjk, < j,k < r? , while 
{ , } represents anticommutation. 

The third constraint on 7^* is that they be positive maps, namely that they 
transform positive matrices into positive matrices. Indeed, as a consequence of the 
duality relation (13), the 7^ are positive maps if and only if their dual maps preserve 
the positivity of states p G S{S); more precisely 

0< Tr(7t[p]x) ^Tr(p7;[X]) >0. 

It has been stressed that positivity preservation is necessary for physical con- 
sistency, but it is not sufficient against couplings with generic ancillas: for that the 
stronger request of complete positivity is compulsory. 

The Hamiltonian contribution to the generator guarantees complete positivity 
since it gives rise to the standard unitary evolution U( that turns out to be auto- 
matically of the simplest Kraus-Stinespring form. Evidently, positivity and complete 
positivity both depend on the properties of the Kossakowski matrix . 

Remark 3.2 Asking for positivity preservation results in quite an intricate alge- 
braic problem; in fact, no general necessary conditions, but only sufficient ones are 
so far available on the coefficient Cij such that they give rise to positivity preserving 
semigroups 7t (for more details, see Refs. [56-59]). On the contrary, the condition 
that the 7t be completely positive is far more stringent. □ 

Theorem 3.2 The semigroup {'yt}t>o consist of completely positive maps if and 
only if the Kossakowski matrix is positive definite. 
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Remarks 3.3 

(1) The proof of Theorem 3.2 was given in Ref.[54] for finite-dimensional sys- 
tems and in Ref.[55] in any dimension, under the assumption of boundedness 
of the generator. Generators with positive Kossakowski matrix arc known as 
Kossakowski-Lindblad generators and the resulting semigroups are known as 
quantum dynamical semigroups. 

(2) Using the duality relation (13), one gets the following master equation on the 
state-space: 

L[p] = -i{B, p\ + ^ cJf]pF, -[f^f] , p} j . (49) 

The generated dual maps 7t preserve the hcrmiticity and the trace of p; if 
the Kossakowski matrix \Cij\ is positive definite then 7( ® id^ preserves the 
positivity of all states of S and of the compound systems 5 -|- ^4 for any choice 
of finite-dimcnsional^ancillas A. 

(3) By setting K = Yli,j^i ^ij ^iFJ , together with the Hamiltonian the anticom- 
mutator in (49) can be incorporated in a pseudo-commutator: 

-i[H, P\ - III Cij{F,F] , p} = -i^H -Ik^P + ip[H+'-K^ . 

The latter is the typical phenomenological expression for the generator of the 
time-evolution of a decaying system, where K describes loss of probability that 
is irreversibly transferred from the system S to the decay products. 

(4) Beside the pseudo commutator, the remaining contribution to dissipation comes 
from a term that, in the case of completely positive 7(, can be put in Kraus- 
Stinespring form as done in (30). Such a term corresponds to what in classical 
Brownian motion is the diffusive effect of white- noise. It is of the form of 
the wave-packet reduction mechanism (16) and it is interpreted as quantum 
noise. □ 



The consequences of the previous two theorems are best exposed in the two- 
dimensional setting as showed in the next example. 

Example 3.1 Let n = 2 as in Examples 2.2; in such a case, by choosing the 
orthonormal basis of Pauli matrices Fj = j = 1, 2, 3, the dissipative contri- 

bution to the semigroup generator reads 



where a factor 1/2 has been absorbed into the coefficients Cy. In the two- 
dimensional case, it proves convenient to adopt a vector-like representation; we 
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write density matrices in the form 



I2 + p - 



1 
Pi 



P3 

- ip2 



Pi 
1 



IP2 
P3 



, 0<Det[p] = -(l-^pn ,(51) 



where a ~ (cri, (T2, cts) and pis a vector in R^, of unit length if and only if p is a pure 
state; it is usually referred to as the coherence^ or Bloch* vector. By representing 
p as a 4- vector \p) = (1, pi, p2, Ps), any linear operation p i-^- A[p] corresponds to 
a 4 X 4 matrix £ — [£^1/] acting on |p). Thus, the evolution equation (45) can be 
recast in a Schrodinger-like form (the —2 in front is for sake of convenience) 



dt\pt)= -2{H + V)\pt) , 



(52) 



where the 4x4 matrices H and V correspond to the commutator Ljj and to the 
dissipative contribution D. 

It is no restriction to take the Hamiltonian of the form H = u) 
and LJ = (wi,a;2,W3) e M^; using the algebraic relations [ai , aj] = 
i,j = 1, 2, 3, it follows that 



• (T, with 
2* El 





ij k •) 



/o 








\ 








UJ3 


—UJ2 












Vo 


U}2 




/ 



(53) 



This is the typical anti-symmetric form of the action of commutators —i[H , •] when 
represented as a matrix acting on |p). 

Concerning the dissipative matrix V, the requests of trace and hermiticity preser- 
vation impose T>oj = 0, j = 1,2,3, and V^i, e M. By splitting V into the sum of 
a symmetric and antisymmetric matrix, after incorporating the latter into H, one 
remains with 



V 



u 

V 




b 
a 



c 
P 

7/ 



(54) 



where the nine real parameters depend on the phenomenology of the system- 
environment interaction. By exponentiation, one gets a semigroup of 4 x 4 matrices 
Qt = e~2*('^+^); its correspondence with the semigroup {7f}t>o on the state-space 
S{S) is given by 



p p{t) = 7t[p] = ^ p^(i) cr^ 



where p^{t) are the components of the 4-vector \pt). 

Since the trace is preserved at all times, checking positivity preservation amounts 
to checking whether Det[p(t)] > for alH > and for all initial p. The contributions 
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of the anti-symmetric Ti cancel out, thus only the dissipative term remains and the 
time-derivative of the determinant reads 

■ 3 3 " 

^ Vijpipj + ^Vjopj . (55) 

Li,j=l j=l 



dDet[p(t)] 



At 



t=o 



Let p be a pure state: P{n) = (I2 + n- a) /2, with n a unit vector, ||n|| = 1; then 
Det [P(n)] = 0. For positivity to be preserved, it is necessary that 



b[P{n)] 



3 



> . 



Sending n \ 



—fi, the same argument for the pure state P{—n) yields 
3 3 



D[P{-n)] = 2 ( ^ VijUiTij - Y^Vjonj) > . 



Summing the previous two inequalities and varying n in the unit sphere, it turns 
out that positivity is preserved only if 

(a b c\ 
b a ^ > . (56) 
c p 

The above condition on V'-'^^ is only necessary to positivity preservation as 
D[P{n)] < can follow because of the presence of the additional contribution 
J2^=i^joPj- However, it becomes also sufficient when we ask that Qt does not 
decrease the von Neumann entropy of any initial state, as this is equivalent to 
u = V = w = in v. Indeed, the totally depolarized state I2/2, i.e. (1,0,0,0) in 
vectorial notation, has maximal von Neumann entropy log 2; if we impose that 
Qt does not decrease its entropy, it can only stay constant, which means that 
I2/2 is a stationary state. In general, the vectorial expression of stationarity reads 
{H + 'D)\p) = 0, which in the case at hands implies u = v = w = 0. 

In the following, we shall restrict to entropy-increasing semigroups; then, in 
terms of the entries Cy , the matrix V in (54) reads 

/O \ 

C22 + C33 — C12 —Ci3 

— C12 Cii + C33 — C23 

\0 — Ci3 —(^23 C'li + C22 J 

Thus, the positivity of [Cy], which, according to the previous discussion, is necessary 
and sufficient for the complete positivity of 7t, results in the inequalities 

2R = a + j-a>0, RS > b'^ , 

2S = a + 'y-a>0, RT>c^, 

2T = o-|-Q!-7>0, ST>0^, 

RST>2bc(3 + Rp'^ + Sc^ + Tb'^ . (58) 



V = 



(57) 



26 Benatti, Floreanini 



These constraints are much stronger than those coming from positivity alone, that 
is from D^^^ > 0, which yields 

a > ( aa>b'^ 

a>0 , < a7 > c2 , DetD^^^ > . (59) 
7 > i a7 > 

□ 



The fact that the conditions for complete positivity arc stronger than those for 
positivity have an important physical consequence in that the decay-times related to 
dissipative completely positive semigroups must obey a definite hierarchy as showed 
in the following example. 



Example 3.2 A typical relaxation behaviour induced by a heat bath on a two 
level system is determined by the following evolution equations for the entries pij 
of its density matrix: 



dpii 

dt 

dp22 

dt 
dpi2 

dt 
dp2i 

dt 



= -P Pii 



q P22 



= -q P22 + P pii 
= -{ico + r)pi2 , 

= (iuj - r)p2i , 



(60) 
(61) 
(62) 

(63) 



where p, q, r and iv are positive constants, the latter representing the Hamiltonian 
contribution. The above time-evolution equations have been written for the stan- 
dard matrix representation p = { ) ; going to the vectorial representation 

V P21 P22 / 

introduced in the previous example one gets 



dt\pt) = - 
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p + 


J 







(64) 



The choice p = q induces entropy-increase, so that positivity-preservation is equiva- 
lent to the positivity of V^'^'' which only requires r,p> 0. On the contrary, from (58), 
complete positivity asks for r > p, i.e. the so-called phase-relaxation I/T2 = r must 
be larger or equal than l/(2Ti), where 1/Ti is the so-called population-relaxation. □ 



Remark 3.4 The hierarchy between characteristic decay times induced by the 
request of complete positivity is often refused or challenged in the physical litera- 
ture {e.g. see Refs.[39, 60-66]). Indeed, as already stressed, it arises because of the 
possibility of an uncontrollable statistical coupling of the experimentally accessible 
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system with an inert entity, the anciUa. Such abstract eventuahty seems to be too 
weak to constrain the physical properties of dissipative quantum systems. We shall 
later show how the true role of complete positivity is far from being an abstract 
mathematical non-sense, but has to do with the possibility of entangled bi-partite 
states evolving dissipatively in a same environment. □ 



3.2. Master Equation 

In this Section, we describe how from the global time-evolution equation 

dtPS+E{t) =l.s+E[pS+E{t)\ , (65) 

by performing certain Markov approximations, one can obtain master equations as 
in (45).^'^'"^ We shall be particularly interested in comparing them with the general 
form in (49). 

The first step consists in defining a projection P : S{S + E) ^ S{S + E) on 
the state-space oi S + E which, by partial tracing, extracts the state of S: 

P[ps+e] = TiEips+E) O Pe , (66) 

where is a chosen environment reference state. It follows that P"^ = P o P = P 
and, with Q such that P + Q = 1s+e, 

P[ps O Pe] =Ps<S)pE , Q\ps O Pb] = . (67) 

The second step uses P and Q to split the evolution equation into the two 
coupled equations 

dtP[ps+E{t)] = ^S+E[P[PS+E{t)]] + l.'s^lQlps+Eit)]] , (68) 
dtQipsMt)] = Lg^^[P[ps+B(i)]] + L^%[Q[ps+E{t)]] , (69) 

where L^^^ = P o Ls+e o P, = Po Ls+e o Q, L^^^ = Q o L5+B o Q and 

Lf = Q o hs+E o P. 

Formal integration of the second equation yields: 

Q[ps+E{t)]=c'^s?.[Q[p,^^(0)]] + rdse(*-«)^??-oLf^^[P[ps+ij(s)]] .(70) 

Jo 

We shall now assume an uncorrelated initial state of the tensor product form ps®Pe- 
Thus, the first term on the right hand side vanishes because of (67); once inserted 
into (68), the solution (70) leads to the following evolution equation 

dtPs{t)®PE=l^'slE[ps{t)®PE] + /'dsLf^goe(*-^)'^s?-oLf^^[ps(s)®pB] .(71) 

Jo 
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The tensorized pe on the left hand side also appears on the right hand side because 
of (67)-(69); it can thus be traced away yielding a master equation for ps{t) = 
Tr E[ps+E{t)]. 

At this point, one assumes the reference state to be an equilibrium state of 
the dynamics of the environment alone, 

LeIpe] = -i[He , pe] = ■ (72) 

Such a condition is physically plausible because typical environments are mod- 
elled as very large heat baths or thermostats whose equilibrium states are not per- 
turbed by the interaction with S. As a consequence, (I5 ^•e)P = 0; similarly, 
P(ls ® 1L_e) = follows from the cyclicity of the trace operation and expresses the 
conservation of probability in the environment. 

Further, in (17), one assumes an interaction Hamiltonian of the form 

H' = J2ya^Ba, , (73) 

a 

where and S„ are hermitian operators acting on the Hilbert spaces C" of S, 

respectively Sj of E. It proves convenient to define centered environment operators 
Ba = Boi — TiEiBa); in this way one gets a new Hamiltonian for S and a new 
interaction term 

H^ = Hs + XY^^c. TTEiB^) , H' = J2Va®B^ , Tt[pe B^] = . (74) 

a a 

^ V ' 



Remark 3.5 The redefinition of Hs amounts to a Lamb shift of the energy levels 
due to a mean-field, first order (in A) approximation of the interaction with E. □ 

Denoting by L;^ (— i times) the commutator with respect to Hg and by L' (— z times) 
the commutator with respect to the new interaction term H' , the total Liouvillian 
in (65) becomes 

L = hs+E = «) lis + Is «> Lb + AL' . 
Moreover, by means of (72) and of the last relation in (74) one gets 

Po (L^(g)l£;) = (L^0l£;)oP , PoL'oP = 0. (75) 

Using the latter relations and computing the partial trace Tr^; with respect to the 
orthonormal basis of eigenvectors oi pe, further yield 

^s+E[ps{t) (g, pe] = (L^ ® lEf^[ps{t) pe] = l^siPsit)] <^ PE , (76) 

l^^.%[ps+Eit)] = A (l.Y'^iPs+Eit)] = ATrB(L' o Q[ps+Eit)]) ® Pe , (77) 

i^1+E[Ps{t) ® Pb] = A (Vflpsit) ® pe] = Xh'[ps{t) ® pe] , (78) 
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where the last equality follows from 

QL'[ps ® pe] = L'[ps ® pe] ■ (79) 

Inserting these expressions into (71) and getting rid of the appended factor pE, one 
finally gets the master equation 

dtps{t)=h^s[ps{t)] + ^*dsTrB(L'oe(*-^)"^^^ o L' [ps(s) ® Pe]) • (80) 

The second term in the right hand side of (80) compactly comprises the effects 
due to the environment; these are described by two nested commutators with re- 
spect to the interaction Hamiltonian H', with the full time-evolution present in 
exp((t — s)L'3'3), and through a final overall trace over the environment degrees of 
freedom. Finding the concrete form of these effects would require the solution of 
the time-evolution of the total system S + E. 

In order to proceed to a more manageable equation and to the elimination of 
all memory effects, one first notices that Tr^ makes the dissipative term ultimately 
depend on two-point correlation functions of the environment operators Ba ■ Mem- 
ory effects are thus expected to become negligible on a time-scale much longer than 
the typical decay time te of correlations in the environment. 

That indeed this is the case can be illustrated with the help of the following 
simple, exactly solvable model of system-environment coupling.^'' 



Example 3.3 In the notation of Example 3.1, consider a two- level system im- 
mersed in a bosonic thermal bath at inverse temperature (3 = T~^. The Hamiltonian 
for the total system contains three contributions as in (17): 

H = ^0-3 (g) 1e -I- Is '^^k alak + (73 ^ (Xk al + Afe ak) , (81) 




where, for simplicity, the bath modes are assumed to be one-dimensional. The first 
two pieces, giving together the free Hamiltonian Hq, drive the uncoupled motion: 
Q represents the splitting between the two system levels, while Uk the energy of the 
bath modes, labelled by the discrete index k (the bath is for the moment assumed 
to be confined in a finite "box"); the corresponding creation, a]., and annihilation, 
Uk, operators are taken to satisfy the standard oscillator algebra: [a|.,a/c'] = 6k,k', 

[«L4'] = [afe,afe'] = 0. 

The coupling between the subsystem and the bath degrees of freedom is de- 
scribed by the third contribution in (81), with Afc, Xk playing the role of mode- 
dependent coupling constants. Note that the interaction Hamiltonian H' is diago- 
nal in the subsystem Hilbert space: it is this property that allows deriving an exact 
analj^ic expression for the total evolution operator Ut = e~**^. 

In order to get rid of the free motion, it is convenient to work in the interaction 
picture. Then, after conjugation with the free evolution operator U^^^ = e"'*^". 
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Ut can be obtained through the time-ordered exponentiation of the foUowing time- 
dependent interaction Hamiltonian: 



+ Afe e 



(82) 



k k 

Since the operators aI. = as ^ a]., Ak = as ^ Uk obey the same algebra as a|., au, 
this exponentiation can be computed in closed form using standard Lie algebraic 
manipulations.*'^^ As a result, the time-evolution operator in the interaction picture 
reads (T represents time ordering): 

U{to,t) = Te-'-^l ''^"'^^'> = nkTc-' ll '^-"'"(-^ 



where 



i^{t - fo) = ^ (^kit - to) - sina;fe(t - to) 

Kito,t) = J2 [Mto, t) 4 - Ik{tQ,t) G 



(83) 



(84) 



and 



fk{to,t) 



drAfc e«^'=^ = — ( c*'^'=*« 



(85) 



At the initial time to, subsystem and environment are assumed to be uncorre- 
lated, with the bath in a thermal equilibrium state described by the density matrix 
pE = e~''^®/Tr[e~''^^]. The evolution in time of the reduced density matrix p{t) 
pertaining to the subsystem is then given by: 



p{to) p{t) = Tte \u{to, t) p{to) ^ PE U^ito, t) 



(86) 



A convenient basis in the two-dimensional subsystem Hilbert space is provided by 

the eigenvectors of Hs, o^li) = (^l)*^^N)j * = 0, 1. Using (83), from (86) one then 
deduces the evolution of the corresponding matrix elements Pij{t) = {i\p{t) \ 



Pooit) = Poo (to 
Pio{t) = Tte 



Pll{t) = Pll{to) , 

Pio(to) =e-^(*°'*)pio(to) , 



with 



T{to, i) = 4 V M [1 _ cosa;,(t - to)] coth ( ^) 

u ^k \ ^ J 



(87) 



(88) 



Only the off-diagonal matrix elements are affected by the interaction with the bath: 
in this case, the environment is responsible for the loss of quantum coherence, 
without affecting the population of the two levels; this is a consequence of the form 
of interaction Hamiltonian H' , which is diagonal in the chosen system basis. 
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In order to analyze in detail this decoherence phenomenon, it is useful to pass 
to the thermodynamical limit of the bath system: in this way, the mode label k 
becomes a continuous variable. The damping function (88) can then be expressed 



as an integral over the frequency variable (^j, (l/27r) duj) 



T{to,t) = T{t-to) 



2A2 



1 — cosu>{t — to) 



coth 



(89) 



where a typical dependence of the coupling constants on frequency has been as- 
sumed, together with a massless dispersion relation. The introduction of the damp- 
ing exponential term is necessary to make the integral converge at infinity: physi- 
cally, it is justified by assuming the presence of a characteristic scale 1/e below which 
the system-bath coupling decreases rapidly, as found in many realistic models. 

The function in (89) can be conveniently split as V(f) = To{t) + r^(t) into a 
vacuum, ro(t), and thermal, I'pit), contributions, that can be separately evaluated. 
The vacuum piece, 



ro(i) 



_ r 

~ 7!" io 



do; 



1 — cos Lot 



cotM^I-i 



TT 



l + {t/e) ,(90) 



is temperature independent and describes how the bath vacuum fluctuations aff'ect 

the quantum coherence of the two-level system. On the other hand, the action of 
the bath thermal fluctuations on the subsystem dynamics is determined by^® 



4A2 



ds 



do; sin ut 



--—In 



\r{l + e/p + it/f3)\ 
r(l + e//3) 



(91) 



The existence in the problem of different time scales allows to identify, through 
the explicit results (90) and (91), three time regimes in the two-level system dy- 
namics: 

(1) short times, t <C e, for which T{t) ~ (A^/7r)(t/e)^; in this case the bath fluctu- 
ations hardly affect the subsystem free evolution; 

(2) intermediate times, s < t < P, where r{t) ~ (2A^/7r) ln(t/e); here the main 
cause of decoherence is due to the vacuum fluctuations; 

(3) long time regime t 3> 0, where r{t) ^ (2\'^//3)t.^ This is the M^arkovian regime, 
which holds for times much longer than the characteristic decay time of the bath 
correlations, which is proportional to the inverse temperature, te — /?. The loss 
of coherence follows an exponential decreasing law, with a lifetime T2 = /3/2A^; 
this corresponds exactly to the phase relaxation parameter introduced in Exam- 
ple 3.2. As already mentioned, the population relaxation parameter Ti is here 
vanishing due to the particularly simple choice of subsystem-bath coupling. 

□ 



''This result follows from the following asymptotic behaviour of the modulus of r{x + iy) as 
\y\ 00, for X and 2/ real: |r(a; + iy)\ ^ V^n |2/|^-i/2 e-^\y\/2 69 
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3.3. Markovian Approximations 

Obtaining a physically consistent Markovian approximation for the master equation 
in (80) is notoriously tricky.^'^'^'™ Indeed, with respect to the previous Example, 
the situation is complicated by the dynamics of the system S alone, generated by 
the Hamiltonian Hs , which in general does not commute with the interaction with 
the bath. This introduces in the problem another timescale rs, and the need of 
a hierarchy condition te <^ rg, that allow a clear separation between subsystem 
and environment. Following Rcf.[71], wc shall now present the most commonly used 
Markovian limits of (80), focusing in particular on their physical consistency, espe- 
cially in relation to the property of complete positivity. 



3.3.1. Second Order Approximation 

One of the most simple assumptions allowing an immediate identification of sub- 
system and environment is the hypothesis of weak coupling between the two. In this 
case, the ratio te/ts is small because rg — + oo, while te remains finite, and further 
A -C 1. The master equation (80) can be more explicitly rewritten as 



dtPs{t) = -i Hs + XH^i\ ps{t) 



H' ,ps{s)®pE 



(92) 



making clear that dissipation becomes relevant on a slow time-scale of order A~^. 

A first Markovian approximation to (92) can be euristically obtained as follows; 
integration of (92) yields: 



ps(0=e*^-[ps(0)] - A' 



f\v f 
Jo Jo 



X Tte ( [h' , ef"-") "^s^E ^H' , PS {u) ^ Pe]]) ■ (93) 

Then, by changing the integration order and by introducing the new integration 
variable w = v — u, (93) becomes: 



ps{t)=e''^s[psm - X' f due^'-^'"^^ x 

Jo 



dwe-"''^sTr 



H' 



H' , psiu)® pe 



(94) 



In view of the assumed smallness of the coupling constant A, the dissipative 
character of the dynamics emerges on the slow time scale t — }?t. A first instance 
of Markovian approximation consists in 1) substituting r/A^ for t in the curly 
bracket and 2) in letting A — > in the upper limit of the integration in w and in 
the generators and 'Ls+e- As 'hs,E commute with Q and because of (79), this 
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finally gets 

ps{t) = e'^^Hpsm - >? tdu e(*-")"^s x 

Jo 



which is the formal solution to 



dtpsit) 



jji gtu(Lffg+LE) 



H' , ps{u) (g) PE 



where 



L [ps] = - J^ds TvE {e-'^"s , e« (^s+l^) ^ ^ 



(95) 
(96) 
(97) 



Master equations of the form (96) with the dissipative contribution as in (97) are 
known as Redfield equations}^ Given the form (74) of the S — E interaction, we set 

Vo,{t) = e-*'^^[V;] = e'*^^ 14 e-'*^« , B„(t) = e-*'^^[B„] = e'*^^ e"'*^^ .(98) 

Further, we introduce the environment two-point correlation functions 



Gap{s) = Tr pEBa{s)Bp 



Tr 



pE BaBfj{-s) 



(99) 



The dissipative term Di[ps(t)] exhibits the following structure of the environment 
induced second order effects: 



©i[ps(i)] 



V dJ G^0{s) k(s) , V0ps{t) 



+G0a(-s) Ps{t)V0 , ya(s) \ . (100) 



In order to compare them with the Kossakowski-Lindblad generator (49), it is 
convenient to introduce the operators 



C| = V/ dsG^p{s)Vc,{s) 

n, Jo 



(101) 



and rewrite 



E'i[Ps(i)] =l^{Bf ps(i)C| + (C|)^ ps{t)Bl 


-qf B^f ps{t) - ps{t) Bfj {Cfj^} . (102) 

It can be verified that trace and hcrmiticity arc preserved, namely that 
TT(Bi[ps{t)]^ = and that (Bi[ps{t)]j = ]D)i[/9s(t)]. Despite a certain similar- 
ity to (49), in order to check whether the semigroup generated by the Redfield 
equations consists of completely positive maps one has to extract from (102) the 
corresponding Kossakowski matrix and ensure its positivity. Instead of doing this. 
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in the following example we shall show that, in general, the semigroups generated 
by the Redfield equations do not even preserve positivity. 

Up to this point we have considered any equilibrium environment state ps', 
however, typical environments are heat baths and pE thermal states at inverse 
temperature /? = T~^. They satisfy the so-called Kubo-Martin-Schwinger (KMS) 
conditions'"^ 



Ga0{t)=G(3^{-t-i(3) . 



(103) 



These equalities express the analiticity properties of thermal correlation functions 

with respect to time; they arc easily derived when He has discrete spectrum, but 
survive the thcrmodynamical limit and then hold for truly infinite environments. 

Example 3.4 As in the previous Example, consider a two-level system immersed 
into a free Bose-gas in equilibrium at temperature T = /3~^. This time, however, we 
choose an interaction Hamiltonian which does not commute with the system one: 



Hs+E = ^0-3 OIb + I5 (8) ^ LOkoiak +Xai<^B ; 



2 

Hs 



(104) 



H' 



He 



here, B = is any observable of E with Ty[pe B] = 0, where, as in Example 3.3, 
the environment reference state is given by 

q-I3He 

Using the fact that 

C7i(s) =e^""'^^/2£7ie-^""'^^/2 ^cos(sn)(Ti - sin(sn)o-2 , 

the resulting master equation derived by means of the same Markov approximation 
used for (100) reads, with G{s) = Ti[ pE B{s)B], 



dtPsit) = i^ [(73 , Ps{t)] + A' Bi[ps{t)] , 



(106) 



where 



i-[ps{t)]= J ds|G(s)(^cos(sn) (Ti,aips{t) - sin(sf2) (T2 , cti ps{t) ^ 



G{—s)(^cos{sQ) ps{t)o-i,ai — sin{sQ) ps{t)cri,(J2 



(107) 



In order to discuss the implication of this Markov approximation, it proves conve- 
nient to adopt the vector representation (52). Therefore, we insert the expansion 
of ps{t) along the Pauli matrices as in (51); this yields the following expression for 
the matrix 7i and V: 



n = 






\o 



0\ 

n/2 

-0/2 

0/ 



/O 



\d 




b 
a 







aj 



(108) 
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where 



a = ^ ds cos(sf2) {g{s) + G(-s)) 



b 

d 



\2 r+oo ^2 

2 7-, 



+ 00 



dse 



G{s) 







iX 

\2 /■+00 



ds sin(sO)(G(s)+G(-s)) , 

ds sin(sf7)(G(s)-G(-s)) 



\2 p+oc \2 /•+00 

— / dse'^*G(s) - — / dse-*"*G(s 

2 J-oo 2 7_oo 



(109) 
(110) 

(111) 



Because of b, the 3x3 matrix I)\ 

\0 a 

splits into a symmetric and anti-symmetric part: 



6 0' 

a 1 is not hermitian. However, it 







0\ 

6/2 

-6/2 

VO 0/ 



/O 0\ 

6/2 

6/2 a 

\d aj 

V 



(112) 



the contribution Ti*^^^ can be absorbed into 7i as a second order correction to the free 
Hamiltonian. Now the corresponding matrix I?'^-* is hermitian, but, again because 
of 6, apparently non-positive. Suppose 6 = 0, in order to check that the generated 
semigroup is positivity preserving, according to (59), it is necessary that a > 0. 
This is indeed so: in fact, from the KMS conditions (103) apphed to the thermal 
state pb, it follows that 



/+00 
ds e'"" G{s) , a - d = e->^^{a + d) 
-OO 



(113) 



If we thus show that a + d> 0, then a > 0; this is achieved by explicitly computing 
the two-point correlation function within the simplifying framework where is 
effectively a density matrix, with orthonormal basis HeIsj) = Sj\ej). With Zp the 
normalization oi pE, 



a + d = Zfj 



-t-oo 



dse'^^Trf e-0HE^isHE Be-'^^^B 



/ + OC 

= 27r^^ ^e-^^^- |(e,-|S|£fe)|2 5(0 + Sj - efc) > . 



(114) 
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We are now in position to study the consequences of the presence of b and d. We 
shall use the same argument developed in Example 3.1; at time t = 0, (55) gives 

^Det(p) = 26pip2 + apl + dps + ap% . 

One checks that 

1 /4a2 _ (p h /4a2 _ rf2 ^ 

fulfil Pi+ P2 + p\ = ^ and thus give a pure state p = P (Det(P) = 0). Moreover, 
unless 6 = = 0, 

— Det P = - o , ,0/ < • 
At 4(a2 + }/■) 

Because of (103), even if & = 0, = is possible only at infinite temperature, that 
is when /3 = 0. It thus follows, that at finite temperature, the Markov approxima- 
tion discussed in this Section does not preserve positivity, since a state like P is 
immediately turned into an operator with negative determinant, thus signalling the 
emergence of negative eigenvalues in its spectrum. 

Despite this physical inconsistency, the semigroup 74 generated by Ti+V in (108) 
is nevertheless often used in phenomenological applications on the basis of its good 
asymptotic behaviour. In fact, it has a unique equilibrium state peq which is easily 
obtained from the condition (W + 'D)|peg) = 0, yielding pi= p2 = Q and pz = —d/a. 
Using (113), this gives: 

1 /e-^"/2 

Therefore, as expected, the bath drives the embedded open system at equilibrium 
at its own temperature. Nevertheless, as we shall see next, there is no need to use an 
unphysical treatment in order to get a Markovian evolution leading to asymptotic 
relaxation to equilibrium. □ 

Remark 3.6 The origin of the lack of positivity-preservation lies in the way the 
free time-evolution generated by Hs interferes with the various truncations. Indeed, 
if ri = 0, so that there is no free evolution of S, then d = b = 0. Observe, however, 
that in this case, according to (58), the dynamical maps 7t turn out to be completely 
positive and not only positive. In the next Section, we shall sec a particular way to 
deal with the free evolution that automatically leads to a semigroup of completely 
positive maps. □ 

3.3.2. Weak Coupling Limit 

In order to obtain a physically consistent Markovian approximation of (92) in the 
limit of weak couplings, a more careful derivation is needed. The formal integration 
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of equation (96), 

Jo 



(116) 



can be rewritten as 



e-^'^s ps{t) = ps{Q) + y j ds (^e-"^s oDioe'^^sj oe-"^s[pg(s)] . (117) 
Upon passing to the slow time r = A^t, it reads 

+ j\s |^e-(«/^')^^ oDi oe(^A')^^^ oe-(^A')^^[ps(s/A')] • (118) 
Using (98) and (100) one has 

pOO 



+G/3a(— ") 



a,/3 ' 



y„(u)e-('*/^^^s , e('*/^^^s e-('*/^^^s 



PS 



p^e^-/^')^^ l/^e-(-/^')^^ , e(-/^')^^ 14(«)e- 



(119) 



For sake of simpUcity, we assume now the effective Hamiltonian Hg to have discrete 
spectrum = YZ=\ Pa - By inserting YTa=\ = into (119), we get 



_n_ /.oo 



X< G«/3(u) 



a,/3 a,6,c,(i=l * 



+G)3a{-u) 



PsP^VpP,^ , P,-^ya(«)Pc 



(120) 



When A ^ 0, only the exponentials with e^ — e^ = expected to effec- 

tively contribute, since otherwise the the term within parenthesis in (118) rapidly 
oscillates, giving vanishingly small contributions; this approximation is known as 
rotating-wave approximation. Also, under suitable assumptions, when A — > 0, the 
eigenvalues and the cigenprojections P^ are expected to go into eigenvalues 
and eigenprojectors Pa of Hg- 

In more mathematically precise terms, one can show that the oscillating term 
can be substituted by the following ergodic average 

1 /■+'^ 

D= lim — / ds e-'^^-^s o Di o e^"^^ . (121) 
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Further, when A ^ 0, D commutes with e ^^/■^^'>^s^ so that, one finally gets the 
Markov approximation 

ps{t) = e*^^ [psm + /* ds e(*-^)"^^ o B[ps{s)] , (122) 

Jo 

which solves the master equation 

dtpsit) = l^s[ps{t)] + X'mpsit)] . (123) 

Before discussing in details the form of the dissipative term ID), we reconsider Ex- 
ample 3.4. 



ExEunple 3.5 In the vectorial representation of Example 3.1, the free time- 
evolution amounts to the 4x4 matrix 





/I 





Vo 



cos(,sri) 

sin(sr2) 






— sin(sr2) 
cos(sn) 







1/ 



which yields the ergodic average 



lim 



2TJ_t 



( 



\ 






\d 




b 
a 




0\ 




aj 



J 






\d 




a/2 
-6/2 










6/2 





a/2 








a) 



(124) 



From the averaged matrix we extract the antisymmetric part 






Vo 






-6/2 






6/2 








0/ 



(125) 



which gives a second order correction to H in (109), and the dissipative term 









a/2 







a/2 







aJ 



(126) 



A part from the entry d, we see that the 3x3 matrix Dg^^ satisfies the conditions (58) ; 
also, the presence of d docs not spoil the complete positivity of the generated map 
7t, for, as we shall show, the ergodic average automatically yields completely pos- 
itive semigroups. We conclude the example by noticing that the equilibrium state, 
determined by {H + W^^) + V2)\peq) = is as in (115). □ 



In order to prove the claim at the end of the preceding Example, we now study in 
detail the structure of (121) and show that it amounts to a generator as in (47)-(48) 
with positive Kossakowski matrix. 
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We set Lo = Ea — £b with £„ running over the spectrum of Hs ; by means of the 
decompositions 

V„{lo)= PaVaPb, V0{-u;)= J2 PcVpPd = V^{uj) , (127) 

ea—Sb=Ul e<j— £o=W 

from (120) and (121), we expUcitly get 



a,/3 u> Jo I 



Vo,iL0),V^{L0)ps 



Further, it follows from standard Fourier analysis that^ 



where 



/ + 00 
-OO 



is a positive matrix,' while 



Sc./3(w) = —V \ dw — — - = s^„(a;) 



. (128) 



(129) 



(130) 



(131) 



and V indicates principal value. In terms of the Fourier transforms of the two-point 
correlation functions of the environment, one then rewrites 



EE Saf3i'^)Va{oj)V^{l^) , pS 
.a,0 ui 



(132) 



a,/3 w 



Da 



We thus arrive at the following explicit form of the generator (123) 



dtps{t) = -i 



Hs + XH^e'^ +\^Hf ,ps{t) 



(134) 



We thus see that, beside the second order dissipative term of the form (49), the 
environment contributes with a first-order mean-field term and also a second order 



'This can be proved by applying to hag{w) considerations similar to the ones employed in (114) 
of Example 3.4 that use the spectral decomposition of the environment time-evolution. 
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term to the redefinition of an effective Hamiltonian of S. Further, the Kossakowski 
matrices [hajdico)] are positive- definite for all w. 

Remarks 3.7 

(1) Formally, the last statement above follows from Bochner's theorem (see for in- 
stance Appendix A. 3 in Ref.[l]) which states that the Fourier transform hapioj) 
of a function of positive type, as the two-point correlation functions Gaisit), 
gives rise to a positive-definite matrix. 

(2) The ergodic average elimination of the rapid oscillations due to the free system 
dynamics results in a semigroup consisting of completely positive maps. 

(3) An equivalent derivation of the physically consistent Markovian master equation 
(134) can be obtained through the so-called stochastic limit-J^''^^ it provides a 
different, more abstract justification for the ergodic average (121). 

(4) A so-called quantum FokkerPlanck master equation has been derived in [76] ; it 
provides a description of quantum Brownian motion with additional corrections 
due to quantum statistics. □ 

The mathematical consistency of the previous manipulations rests on some as- 
sumptions on the behaviour of two-point correlation functions. In particular, as 
proven in Refs.[72, 73], the following conditions turn out to be sufficient. The first 
one is a condition on the norms of the system observables Va which are coupled to 
the environment, that is X]q ll^all < oo. The second is an assumption on the en- 
vironment which is assumed to be quasi-free, namely all its higher time-correlation 
functions can be expressed as sums of two-point correlation functions. The third 
regards the decay of two-point correlation functions which is assumed to be such 
that there exists > for which 



with a independent of a, (3. 

Remark 3.8 The derivation of the master equation (134) is based on the assump- 
tion of a weak coupling between subsystem and external bath; a master equation 
of the same kind can also describe a different situation, in which the interactions 
between subsystem and environment are rare, but not necessarily weak. The typical 
instance is that of an atom, effectively modeled as an n-level system, immersed in 
a rarefied gas, that represents the environment. The collisions of the atom with 
the gas particles are described by the scattering operator, which involves the full 
interaction Hamiltonian; however, these are rare events, since the gas is by hypoth- 
esis very dilute. Using mathematically rigorous techniques, one can then show that 
in the so-called low density limit the atom subdynamics is generated by a master 
equation in Kossakowski-Lindblad form.''^'^'^ □ 




(135) 
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3.3.3. Singular Coupling Limit 

Unlike in the weak coupling limit, in the singular coupling regime, it is the decay 
time of correlations in the environment that becomes small, te — * 0, while the 
typical variation time of ps, ts, stays constant. Concretely, the two-point correlation 
functions appearing in the dissipative contributions to the open dynamics are made 
tend to Dirac deltas:^'^-™'™ 



(136) 



so that the dissipative effects manifest themselves on time-scales of order t. 

In order to implement such a behaviour and derive the corresponding master 
equation, we consider (17) and rescale it as follows:^ 



Hs+E = Hs(^1e + e"^ls (^He + e'^ H' . 



(137) 



Here we assume H' = J2a B^, with centered environment operators Ba, i.e. 
Tr[BapE] = 0. Then, the same steps leading to (93) yield 

ps(t) =e*^"s[ps{0)] - £-2 f\iy Tdue^*-")^^ x 

Jo Jo 

X Tve{[h' ,e'''^'"'[H' , ps{v-u)(^^pe]]) • (138) 
By letting u = e^w and £ — > 0, one formally gets 



Ps{t) 



[psm 



du TVs ( \h' , e"^«^ [h' , ps{v) ® Pe\\) . (139) 



Explicitly, the above expression solves the evolution equation 



dtPs{t) 



Hs,ps{t) + mpsit)] , 



(140) 



where 



mpsit)] 



V /'"dsiG„^(s)fy„, V0ps{t) 



+G0ai-S) Ps{t)Vp , 14 



(141) 



As was done in the weak coupling regime, one writes 



f 

Jo 



^ / \ "-0/3 , • 



so that (141) splits into a second-order Hamiltonian contribution, with Hamiltonian 

a,/3 
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and a true dissipative term of the form: 



^3[ps]=Y.hc,Jv0psVa - l{VaVp,ps} 



(143) 




(144) 



Remarks 3.9 

(1) The singular couphng hmit, as well as the weak coupling limit, leads to a gener- 
ator of the Lindblad form (49) with a positive Kossakowski matrix and thus to 
a dissipative semigroup consisting of completely positive maps. Notice that, un- 
like in the weak coupling limit case, the operators contributing to the dissipative 
term Da are hermitian. 

(2) By going to the slow time r = s'^t, the global Hamiltonian (137) becomes 



This shows that the singular coupling limit amounts to a weak coupling limit 
where the free motion generated by the Hamiltonian Hs is of the same order 
of the dissipative effects. ^"'^ 

(3) Unlike the weak coupling limit, thermal Bose or Fermi heat baths provide 
frameworks suited to the singular coupling limit only if their temperature is 
infinite.^'''*'™ Indeed, the two-point correlation functions may tend to Dirac 
deltas in time only if their Fourier transforms tend to a constant, but this is in 
contradiction with the KMS-conditions (103), unless /3 = and thus T = oo. 

(4) The infinite temperature limit is equivalent to add to the system Hamiltonian 
Hs an external Gaussian stochastic potential V{t)) In such a case, the partial 
trace Tr^ with respect to ps, is substituted by the average with respect to the 
Gaussian noise driving V{t)7^ The resulting generator is as in (144) with a real- 
symmetric Kossakowski matrix [ha^] (for a specific instance, see Example 3.8 
below); consequently, the totally depolarized state r„ = l„/n is an invariant 
equilibrium state and corresponds to the state (115) with /3 = 0, that is at 
infinite temperature. □ 

Example 3.6 We now reexamine Example 3.4 in the light of the singular coupling 
limit; thus, we consider the rcscalcd total Hamiltonian 



Hs+E = e"^ Hs <»1e + ® He + eH' . 



(145) 




(146) 



fc 



^ This approximation is very useful in many phenomenological applications, in pai'ticular for treat- 
ing decoherence effects in the propagation of elementary particles in matter; for details, see 
Refs.[81, 82]. 
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where now a = 2 TZe 



dtG{t) 



(147) 



With respect to (108) and (126), in the vectorial representation, the singular cou- 
pling limit gives the following dissipative matrix 

/O 0\ 


a 
\0 a/ 

This result directly follows from (109)-(111) by setting f2 = and thus neglecting 

the free system time-evolution whence h = d = 0. By the same argument developed 
in Example 3.4, it turns out that a > and the generated semigroup consists of 
completely positive maps. □ 



3.3.4. A Convolution-Free Master Equation 

The Markov approximations presented in the previous three sections all start 
from the convolution master equation (80). There is however a derivation which 
is convolution- free and also leads to the Redfield-type equations as in (96).®" 

Within the approximation of weak coupling between system and environment 
(A ^ 1), one starts from the complete master equation for the compound system 
T = S + E, 

dtprit) = (Lo + AL')[pT(t)] , 

where Lq comprises commutators with respect to system and environment Hamil- 
tonians and L' the commutator with an interaction Hamiltonian as in (74), with 
centered environment observables. One then passes to the interaction representa- 
tion prit) = exp{—tho)[pT{t)] and solves iteratively the resulting time-dependent 
evolution equation 

dtPrit) = XV{t)[pTit)] , L'(t) ^ . 

With initial condition pt(0) = Pt(0), the solution can be expressed as a series 
expansion in the small parameter A: 

oo 

PT(t)=X]^'f^fcWi?T(0)] , 
fe=0 

where 

Nfc(i)[pT(0)] = / dSi / 'dS2-- - / ' 'dSfcL'(si)oL'(s2)---oL'(Sfc)[pT(0)] , 

Jo Jo Jo 
and No(t)[/9T(0)] =pt(0). 

Under the hypothesis of a factorized initial state, Pt{0) = Ps ® Pe, with pE 
an environment equilibrium state, Ljj^ [pe] = 0, using the projection introduced in 

(66), one gets 

oo 

P[pT{t)] = TrE[pT{t)] ®PE = J2^'' TrB(Nfe(i)[ps pe]) ® Pe ■ 
Ps(t) 
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Further, setting 

oo 

Mk{t)[ps] = Tr E (mt)[ps <E) pe]) , = J2x''Mk{t) 



(148) 



k=0 



one writes ps{t) = $t[ps(0)] and, by a formal inversion of one gets 

dtpsit) = ^tlfsm = o ^t'lMt)] ■ 

Due to the Tr^; operation, the operators A4k{t) contain fc-timc correlation functions 
of the environment. At this point one makes the assumption of a thermal environ- 
ment for which only even correlation functions matter and, according to the weak 
coupling assumption, expand to lowest order in A^, $t ° — ^^■M2{t). Return- 
ing to the Schrodinger picture and explicitly inserting the environment two-point 
correlation functions one gets the Rcdficld-type equations 



dtPsit) 



Hs,ps{t) +MPs{t)] 



(149) 



where 



a/3 



^A[ps{t)] dsi G^pis) \Vf,i-s) ps{t) , V, 



+G/3«(-s) Vo,, ps{t)Vp{- 



(150) 



The dissipative term D4 in the above time-evolution equation is rather simi- 
lar to (100) and thus suffers from the same problems in relations to positivity- 
preservation. This is explicitly shown by the following example. 



ExEimple 3.7 Let us consider Example 3.4 in the light of the convolutionless 
approach; in the vectorial representation, it gives the following dissipative matrix 





\d 



0\ 



-6 a 

0a 



(151) 



where a, h and d are as in Example 3.4. As done there, we split the 3x3 matrix 
/ 0\ 

2?4^^ = I —b a j into symmetric and anti-symmetric part. The first term 
\ a/ 

provides a second order correction H^^^ to the free-Hamiltonian term, while the 
second one gives rise to a dissipative term 

/O 0\ 

^ ^ -h/2 

^ -6/2 a 

\rf a/ 
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of the same form as Vi in Example 3.4, see the second piece in (112), with h —b. 
Therefore, unless b = d = 0, there exist initial states that, as four-dimensional 
vectors, are mapped by the semigroup 



Gt = cxp 



into vectors, corresponding to matrices with negative eigenvalues. □ 



Further elaborations are then needed to make physically acceptable the dynamics 
obtained by the convolutionless method. 



Remark 3.10 A cure against lack of complete positivity and even of positivity 
has been put forward which is based on the so-called slippage of initial conditions 
(see, Refs. [63-66]). The argument is that any Markov approximation neglects a 
certain initial span of time, namely a transient, during which no semigroup time- 
evolution is possible due to unavoidable memory effects. The transient dynamics can 
be effectively depicted as a slippage operator that, out of all possible initial states, 
selects those which can not conflict with the Markov time-evolution when it sets in. 
As we shall see in the next Section, the slippage-argument may cure the positivity- 
preserving problem, but misses the point about the complete positivity issue which 
emerges in full only when dealing with bi- and multi-partite open quantum systems 
and with the existence of entangled states. □ 



The convolutionless approach is often advocated in situations where, instead of a 
weak coupling to an environment, the open system is affected by a weak stochastic 

potential. As already observed in Remark 3.9.4, one substitutes the trace 
with respect to the environment degrees of freedom by an average over the noise. 
Experimental contexts where one can simulate external noisy potentials seem to 
be feasible; by modulating the noise one can enhance the non-positivity preserving 
character of the Redficld-typc equations and thus expose the physical inconsistency 
of the Markovian approximation on which they are based. ^^'^^'^^ 

In particular, neutron interferometry, which has proved to be an extremely pow- 
erful tool to investigate gravitational, inertial and phase-shifting effects occurring 
inside the interferometer,^ may be used to investigate the notion of completely 
positive open system dynamics. In order to do that, we consider below the case in 
which neutrons while propagating inside the interferometric apparatus are subjected 
to weak time-dependent, stochastic magnetic fields coupled to their spin degree of 
freedom. In this case, the noisy effects result from the classical stochastic character 
of the fields and not from an interaction with a quantum environment. 



''For an account of relevant results, see e.g. Refs.[84, 85]. 
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Example 3.8 The states of the neutrons inside the interferometer can then be 
described by a 2 x 2 density matrix R{t). Its dynamics is generated by a time- 
dependent Hamiltonian that takes the form 



Hit) = ^a3 + a- B(i) , 



(152) 



where the system Hamiltonian Hs is as in Example 3.3, while B(t) = 
{Bi{t), B2{t), Bsit)) represents a Gaussian stochastic magnetic field. We assume 
B(t) to have zero mean, (B(t)) = 0, and a stationary, real, positive-definite covari- 
ance matrix (t) with entries 

G.,(t -s) = {B,{t)B,{s)) = G*^{t - s) = G,,(s - t) . (153) 

The time-dependent evolution equation for the density matrix R describing the 
spin degree of freedom splits into 

dtR{t) = {U+U)[R{t)] , (154) 



Lo[i?(t)] 



2 ^^3 , R{t) 



U[R{t)] 



B(t) •(?, R{t) 



(155) 



Because of the stochastic field B(t), the solution R{t) of (154) is also stochastic; 
an effective spin density matrix p{t) = {R(t)) is obtained by averaging over the 
noise. At time t = we may suppose spin and noise to decouple so that the initial 
state is p= (i?(0)) = i?(0). 

The effective time-evolution of p{t) can naturally be obtained following the tech- 
niques of the previous Section, where the trace over the environment degrees of free- 
dom is substituted by the average over the classical noise, Tr^ — > ( ). The resulting 
master equation has the form (150): 

dtPsit) = -i ^(73 , p{t) + B4p{t)] , (156) 



where 



(157) 



The time dependent aj{—s) explicitly reads: 



c^i(s) 



(T, e 



-isHs = \ " 



k=l 



where 



(cos Qs — sin fls 
sin Qs cos Cls 


Inserting these expressions in (157), and extracting the symmetric and antisymmet- 
ric contributions, the master equation results of the form: 



dtp{t) = -i\Hs+H^^^ , Pit)] + Biipit)] 



(158) 
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where H^"^^ — "^^ij^k^i ^ijk i^ij ffc, and P4 is as in (50), with 

3 

2 



1 r 1 

- ^ / ds [G,fc(s) (7fcj(-s) - U,k{s) Guj{-s)\ , (159) 

3 -.00 

G^k{s)Ukj{-s) + U,kis)Gkj{-s) . (160) 

We now study the positivity and complete positivity of the generated time- 
evolution in relation to the decay properties of the stochastic magnetic field corre- 
lation matrix Gij{t). We will consider two representative cases. 

• White Noise 

Let the stochastic magnetic field have white-noise correlations Gij {t — s) = 
Gij 5{t — s), where the matrix Gij is time-independent, symmetric and by assump- 
tion positive-definite. Hence, = 0, and Gij = Gij, so that the Kossakowski 
matrix results positive definite and the generated semigroup is completely positive. 
Further, note that, due to the delta-correlations, all terms A4k{t), fc > 3, in the 
expansion (148) identically vanish, so that the evolution equation (158) is in this 
case exact. ''^ 

• Single Component Covariance Matrix 

Consider a stochastic magnetic field along the x-direction, B(t) = {B{t), 0,0), with 

{B{t)B{s)) = e-^\*-'\ . 

Then, 

_ nB^ ' '\ _ B- 

Unless = the matrix Gij is not positive definite, and therefore the generated 
semigroup is surely not completely positive. Actually, adopting the vectorial formu- 
lation of Example 3.1 and setting: 

2B^\ , B^n 



2A 


n 























the relevant 3x3 dissipative matrix (56) reads 

/ -6 0" 
= \-b a 
\ a, 

which is of the form of the previous Example 3.7; hence, even positivity is not 
preserved. 
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Although apparently formal, these results are not purely academic: indeed, the 
entries of p{t) are directly accessible to experiments, where, by modulating a back- 
ground magnetic field, one may get close to the stochastic properties of the two 
situations just analyzed. □ 

3.4. Why Completely Positive Semigroups? 

In the previous sections we have analyzed three Markov approximations that lead to 
reduced dynamics consisting of semigroups of trace and hermiticity preserving linear 
maps 7t, i > 0, acting on the state-space of an open n-level system 5*. As emphasized 
before, any mathematical effective description of actual dissipative dynamics must 
in the first place preserve the positivity of any initial density matrix. Only in this 
case, the time-evolution turns out to be consistent with the interpretation of the 
state-eigenvalues as probabilities. 

In the second place, one has also to care of possible couplings with inert and 
remote ancillas A and therefore to guarantee the positivity-preserving character of 
the semigroup of maps 74 (g) id/i, too; this is only assured if 74 is completely positive. 

The weak coupling limit with the ergodic average prescription and the singu- 
lar coupling limit give semigroups of completely positive maps and therefore fully 
physically consistent. On the other hand, without ergodic average or through the 
convolutionless procedure one is led to Redfield-type equations which generate dy- 
namical maps that are, in general, not even positivity preserving. 

Remark 3.11 One can always construct semigroups of linear maps that are only 
positive and not completely positive, as seen in Example 3.1 (c/. Eq. (59)). The 
real problem is the derivation of such kind of dissipative dynamics from a micro- 
scopic interaction with the environment. One knows that these semigroups would 
not be acceptable in the case of couplings with ancillas: it is this inconsistency that 
essentially forbids their construction from microscopic interactions.''^'''^ Indeed, all 
known rigorous microscopic derivations of an open quantum reduced dynamics lead 
to completely positive semigroups. □ 

It has also been remarked that complete positivity guarantees full physical 
consistency at the price of an order relation among typical life-times, which, in 
the two dimensional case is expressed by the necessary inequalities (58) (see also 
Example 3.2). It has been observed that from a strictly physical point of view, it 
appears doubtful that possible couplings to remote and inert ancillas might have 
concrete consequences on actual experimental contexts. In view of this, complete 
positivity is often refused as a mathematical artifact, with nice structural conse- 
quences as the Kossakowski-Lindblad form of the generators in (49), but without 
solid physical justifications. 

Leaving aside the difficulty of deriving a positivity-preserving semigroup from 
fundamental system-environment interactions, one may argue that positivity preser- 
vation might suffice to guarantee physical consistency. Indeed, it is enough to ensure 
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the probabilistic interpretation of the spectra of evolving states, while the danger 
of possible couplings to uncontrollable ancillary systems is a far too abstract spec- 
ulation to be of any physical importance in actual experimental contexts. 

However, we have also seen that the true meaning of complete positivity can 
mainly be appreciated in relation to quantum entanglement of which it is just 
another facet. Quantum entanglement requires at least two systems, whereas, until 
recently, the literature on open quantum systems primarily has been dealing with 
single systems immersed in heat baths. 

It appears then clear why the fundamental role of complete positivity does not 
emerge in full in usual open quantum system contexts. On the contrary, the whole 
perspective radically changes when one considers two systems in a same environ- 
ment, as, for instance, two atoms in a heat bath at some finite temperature. The 
main point is that two such systems can be prepared in an entangled initial state; 
the usual unitary timc-cvolution is completely positive by its very construction and 
thus gentle to entangled states, while some dissipative time-evolution may not be 
completely positive and thus potentially conflicting with initial entangled states, 
making negative eigenvalues emerge in the course of the time-evolution. 

More concretely, suppose the reduced dynamics of just one system S in the en- 
vironment E to be given in terms of a semigroup of positive maps 7t. Consider then 
two such systems (both denoted by S) embedded into E, without direct interac- 
tion between themselves. In first approximation, the states of the compound system 
S + S will evolve according to the dissipative semigroups of maps Ft = jt Jt- 

We stress the difference between the two contexts: the one just depicted, where 
the maps 7t (8) 7t describe the dynamics of two systems S in the same environment 
E, and the other, where the maps 7j (g) icIa drive a single open quantum system S 
coupled to a remote and inert ancilla A that has nothing to do with it apart for 
possible statistical correlations (entanglement) in the initial state of S + A. 

Suppose now that complete positivity is a mathematical artifact and may thus be 
dispensed with in actual physical contexts, where no ancilla effectively intervenes. 
Since we want that at least the positivity of the states be preserved, the crucial 
question is then the following: 

Which are the conditions on jt such that = 7t (g) 7t be positivity-preserving on 
the state-space of the open system S + SI 

For n-level systems the answer is given by the following Theorem. 

Theorem 3.3 // a semigroup of positive maps 7t acing on S{S) describes the 

dynamics of an n-level system S, then the semigroup Tt = ^t®lt acting on S{S+S) 
is positivity-preserving if and only if jt is completely positive. 

The conclusion which is to be drawn from the previous theorem is the following. 
If 7t is only positivity-preserving but not completely positive, then F^ = 7t (g) can 
not be positivity-preserving. Therefore, there are states of S -\- S whose spectrum 
does not remain positive under Ft and which thus loose their probabilistic interpre- 
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tation in the course of time. Such states are necessarily entangled as Tt preserves 
in any case the positivity of separable states as in (37). In fact, 



and 7t is positivity-preserving by assumption. It is thus the existence of entan- 
gled states that makes the request of complete positivity necessary. As already 
observed, unitary time-evolutions are automatically completely positive, while dis- 
sipative semigroups are not in general and their physical consistency depends on 
the way Markov approximations are implemented. 

Remarks 3.12 

(1) The argument of Theorem 3.3 gives a concrete explanation of why complete 
positivity is necessary by hinging on the necessity that at least the positivity 

of the timc-cvolution of bipartite systems be guaranteed. If Tt = ® It is 
positive, then 7^ must be completely positive and is completely positive, too. 
Indeed, by using the Kraus-Stinespring representation (29), one checks that 
the composition of two completely positive maps is again in Kraus-Stincspring 
form, and thus completely positive. In the case above, Tt — (7* ® id) o (id(8) 7(). 

(2) If one considers semigroups of maps of the form Tf = 7(^'' ® 7^^^ with two 
different single-system semigroups, then Tt can be positivity preserving with- 
out 7^^'^'' being both completely positive. In this case, Tt is not completely 
positive; in fact, its generator has a Kossakowski matrix which consists of the 
orthogonal sum of those of the generators of 7^^^ and 7j^\ and one of them is 
not positive. Then. Ti is not robust against the coupling of 5 + 5 to ancillas. 

(3) An important ingredient in the proof of Theorem 3.3 is time-continuity and the 
resulting existence of a generator; without continuity, there are counterexamples 
to its conclusions. For instance, take the transposition map T„, which can not 
be continuously connected to the identity. It is not completely positive, but it 
is easy to check that T„ ® T„ is nevertheless a positive map. 

(4) With reference to the slippage of initial conditions argument briefly mentioned 
in Remark 3.10, its usage in the bipartite context would lead to argue that 
the transient effects prior to the Markovian regime are such that all entangled 
states dangerous to a semigroup of non-positive Tt = ® It should somehow 
be eliminated by the slippage operator. Such an operator might plausibly be 
constructed, but it would be adapted to the bipartite contc;xts; by going to 
multipartite systems it should then be re-adjusted, thus making the whole pro- 
cedure become ad hoc and lose the general applicability it should possess. □ 

4. Open Qucintum Systems and Entcinglement Generation 

In this Section we shall apply the results and considerations discussed so far to 
the study of a specific instance of open quantum system, that of atoms interacting 
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with a bath made of quantum fields. As we shall see, this model, besides having 
an intrinsic theoretical interest, underlays many phenomenological applications in 
chemistry, quantum optics and condensed matter physics. 

As remarked before, the use of the open system paradigm in modelling specific 
physical situations requires an a priori unambiguous separation between subsystem 
and environment. In the weak coupling regime, this is typically achieved when the 
correlations in the environment decay much faster than the characteristic evolution 
time of the subsystem: in this case, the changes in the evolution of the subsystem 
occur on time-scales that are very long, so large that the details of the internal en- 
vironment dynamics result irrelevant. It is precisely this "coarse grained" evolution 
that is captured by the Markovian master equations studied in Section 3. 

Independent atoms immersed in external quantum fields and weakly coupled to 
them constitute one of the most common and interesting examples in which the 
open system paradigm finds its concrete realization. Indeed, in typical situations 
the differences between the atom internal energy levels (~ Tg^) result much smaller 
than the field correlation decay constants (~ t^^) so that a clear distinction be- 
tween subsystem, the non-interacting atoms, and environment, the external fields, 
is automatically achieved. 

The atoms are usually treated in a non-relativistic approximation, as indepen- 
dent n-level systems, with zero size, while the environment is described by a set of 
quantum fields, often the electromagnetic field, in a given quantum state, typically 
either a temperature state or simply the vacuum state. The interaction of the fields 
with the atoms is taken to be of dipole type, a well justified approximation within 
the weak coupling assumption. Even in this simplified setting, that ignores all in- 
tricacies related to the internal atom structure and to the full coupling with the 
electromagnetic field, the model is of great relevance both theoretically and phe- 
nomenologically: indeed, with suitable adaptations, it is able to capture the main 
features of the dynamics of very different physical systems, like ions in traps, atoms 
in optical cavities and fibers, impurities in phonon fields (many concrete examples 
are presented in Refs.[l-16]). 

Despite this ample range of possible applications and the attention devoted 
to it in the recent literature, no particular care is generally adopted in the various 
derivations of consistent subdynamics. As a result, time-evolutions that are not even 
positive, of the type discussed in Examples 3.3 and 3.8, have often been adopted in 
order to describe the subsystem physical properties: we have already pointed out 
that physical inconsistencies might be the outcome of such ill-defined dynamics. 

As a relevant application of the theory of open quantum systems dynamics, in the 
following we shall explicitly derive and study in detail the reduced time-evolution 
obtained by tracing over the unobserved (infinite) field degrees of freedom. For 
simplicity, we shall restrict our attention to two-level atoms in interaction with a 
collection of independent, free, scalar fields in arbitrary d space-time dimensions, 
assumed to be in a state at temperature T = l/f3. 

The time-evolution of a single atom will be first analyzed. Following the steps 
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outlined in Section 3, the weak coupling limit procedure will be adopted; it leads to a 
mathematically sound and physically consistent expression for the atom subdynam- 
ics in terms of completely positive quantum dynamical semigroups. Not surprisingly, 
the atom is found to be subjected to dissipative and decohering effects that asymp- 
totically drive its state to an equilibrium configuration that is exactly thermal, at 
temperature T. 

The case of a subsystem composed by more than one atom can be similarly 
treated. Again for simplicity we shall limit the discussion to the case of two, inde- 
pendent two-level atoms immersed in the same environment of thermal quantum 
fields. The corresponding master equation that, in the weak coupling limit, describes 
the reduced time-evolution of the atoms can be obtained by generalizing the tech- 
niques used in the single atom case. It again generates a semigroup of completely 
positive maps, that can be studied in detail. In particular, the asymptotic equilib- 
rium state can be computed explicitly, and in general it turns out to be an entangled 
state of the two atoms, even in the case of a totally separable initial state. 

This remarkable conclusion might appear at first sight rather surprising. As 
observed before, the interaction with an environment usually leads to decoherence 
and noise, mixing enhancing phenomena typically going against entanglement. 

Exeimple 4.1 As a drastic instance of mixing-enhancing behaviour, let us consider 
the so-called depolarizing channel affecting a two-level system; it is described by the 
following Kossakowski-Lindblad master equation 

dtPt = cTiptai + a2Pt'J2 + crsptas - 2>pt ■ 

Its solution is pt = (1 — e~^*) l2/2-|-e~'**p, so that any initial state p, even a pure one, 
with zero von Neumann entropy, goes asymptotically into the totally depolarized 
state I2/2 with maximal entropy log 2. □ 

Since the closer to the totally depolarized state, the less entangled a state is, 
one generally expects that when a bipartite system is immersed in an environment, 
quantum correlations that might have been created before by a direct interaction 
between the two subsystems actually disappear. And this could occur even in finite 
time.*^'^^ The degrading of entanglement by an external environment is clearly a 
curse in quantum information, where one tries to maintain quantum correlations 
long enough to complete a quantum computation. Various error correcting strategies 
have been devised precisely to counteract environment decoherence.^^ 

However, an external environment can also provide an indirect interaction be- 
tween otherwise totally decoupled subsystems and therefore a means to correlate 
them.^^'^°'^^'^^ This phenomenon has been first established in exactly solvable mod- 
els: there, correlations between the two subsystems take place during a short time 
transient phase, where the reduced dynamics of the subsystems contains memory 
effects. 

Remarkably, entanglement generation may occur also in the long time, Marko- 
vian regime, through a purely noisy mechanism:^^ it is precisely this situation that 
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is relevant in the analysis of the dynamics of two independent atoms interacting 
with the same vacuum quantum field. 

4.1. Single Atom Systems: Master Equation 

We shall first deal with a single atom in weak interaction with a collection of free 
scalar fields at temperature T, in d space-time dimensions. As mentioned before, we 
are not interested in the details of the atom internal dynamics. We shall therefore 
model it, in a nonrelativistic way, as a simple two-level system, which can be fully 
described in terms of a two-dimensional Hilbert space. 

In absence of any interaction with the external scalar fields, the atom internal 
dynamics will be driven by a 2 x 2 Hamiltonian matrix Hs, that in the chosen basis 
can be taken to assume the most general form (c/. Examples 2.2, 3.1) 

Hs = '^ri-a, (161) 

where ai, i = 1,2,3 are again the Pauli matrices, Ui, i = 1,2,3 are the components 
of a unit vector, while uj represents the gap between the two energy eigenvalues. 

To conform with the standard open system paradigm (sec (17)), the interaction 
of the atom with the external scalar fields is assumed to be weak and described by 
a Hamiltonian H' that is linear in both atom and field variables: 

3 

H' = Y,^i®^i{x) . (162) 

i=l 

As the atom is taken to be an idealized point-particle, without size (see Remark 
4.1.2 below), the interaction is effective only at the atom position , /i = 0, 1, . . . , c?, 
that, without loss of generality, can be chosen to be the origin of the reference frame. 
As a consequence, all coordinates vanish, x** = 0, ^ = 1, 2, . . . , rf, except for the 
time variable, = t. The operators <I> represent the external quantum fields, and 
are taken to be spinless and massless for simplicity. They evolve in time as free 
relativistic fields with an Hamiltonian if$ whose explicit expression need not be 
specified in detail. Further, we shall assume the atom and the fields to be prepared 
at time t = in an uncorrelated state, with the fields in the temperature state pp, 
as in (105), and the atom in a generic initial state p{Q). 

Example 4.2 Although apparently very simplified, the model presented above 
is of relevance in various phenomenological applications. As an instance, let us 
discuss the case of a two-level atom in interaction with the electromagnetic field 
in an optical fiber through the standard dipole coupling.^^ In this case, working in 
three-dimensional space, the interaction Hamiltonian takes the form: 

H' = -d- E{x) , (163) 

where d denotes the operator representing the atom dipole moment, while E is the 
electric field at the atom position x. 
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In absence of the field, the atom dynamics can again be taken to be driven by 
the generic two-level Hamiltonian (161), with eigenstates |— ) and |+), while the 
electromagnetic field is quantized using the Maxwell action: 

5 = -iy"d3fdtF^,F'^% (164) 

where = d^A^ — d^A^ is the field strength expressed in terms of the electro- 
magnetic 4- vector potential A^. Without loss of generality, we shall choose to work 
in the familiar Aq = Q gauge, so that only the vector potential A will be relevant. 

Since the dipolc moment d is a vector operator with odd parity, one finds: 
(— ) = = 0. Furthermore, one can write: 

{+\d\-) = {-\d\+Y = \d\e'^ u , (165) 

where u is the unit vector that defines the orientation of the atom in space, while 
^ is a phase. As a consequence, the representation of the operator d on the two- 
dimensional atom Hilbert space has vanishing diagonal elements, so that one can 
rewrite (163) as: 

H' = -\d\ E{x) [e''V+ + e-'V_j , (166) 

where E{x) = u ■ E{x) and (j± = cri ± ia2- Further, note that the phase 6 can be 
reabsorbed into a redefinition of the Pauli matrices, a± e^*^c7±, without affecting 
their algebraic properties; in other terms, there is no loss of generality in setting 6 
to zero. 

Let us now exploit the fact that the atom is inserted in an optical fiber. One can 
always choose the direction along the fiber to coincide with the z-axia, the x- and 
y-axis representing the corresponding transverse directions. The spatial geometry 

becomes essentially unidimensional and the vector potential A can be taken to 
be independent from the x and y coordinates. Moreover, in such a geometry the 
electromagnetic energy flows along the fiber and thus only the components of A 
transverse to the fiber are non- vanishing; further, of these only A = u ■ A really 
couples to the atom. 

Taking into account all these conditions, the Maxwell action (164) reduces to 



dtdz 



dtAdtA-d,Ad,A 



(167) 



where S = J dxdy is the cross section area of the flber. By further noting that 
E — —dfA and by introducing the dimensionless rescaled field (j){t, z) = A{t, 2:)\/S, 
the interaction Hamiltonian (166) can be rewritten as: 

H' = \ dt^{0, z) [a+ + a-] , (168) 

where A = plays the role of coupling constant. It is precisely of the form 

(162), with $1 = dfcf) and $2 = ^3 = 0, where ^ is a l-|-l-dimensional free, massless 
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scalar field; its dynamics is generated by the Hamiltonian i?^ that corresponds to 
the action (167): 



(169) 



The scattering theory for atoms in an optical fiber described by interactions of 
the form (168) has been discussed in Refs.[95, 96] by means of a formal mapping to 
the anisotropic Kondo model. Here instead we are interested in analyzing the full 
reduced dynamics of the atom and not only its asymptotic properties. We shall do 
it by turning to the more general interacting Hamiltonian in (162). □ 

The total Hamiltonian H describing the complete system, the two-level atom 
together with the external fields can thus be written as 

H = Hs'»U + ls'»H^ + >^H' . (170) 

It is precisely of the form discussed in the previous Section, with centered environ- 
mental operators, Tr[p^$j] = 0. It generates the time-evolution of the total density 
matrix ptot , via the Hamiltonian equation 

= I^'hM)] , (171) 

starting at t = from the initial configuration: /Otot(O) = p{0) ® pp- As ex- 
plained before, we want to follow the dynamics of the reduced density matrix 
p{t) = Tr$ [ptot (*)] on time-scales that are long with respect to the decay time 
of correlations in the environment. This is can be obtained by suitably rcscaling the 
time variable, t t/}? and then taking the limit A ^ 0, following the mathemati- 
cally precise procedure of the weak coupling limit discussed in Section 3.3.2. 

The resulting master equation for p{t) takes a Kossakowski-Lindblad form (49) 

^ = -*[ffcff,p(t)] +W)] , (172) 

with D as in (50). The effective Hamiltonian TJeff and the coefficients of the 3x3 
Kossakowski matrix that appear in (50) are determined by the Fourier transform 
of the field correlations: 

/+00 
dte'f*($,(t)$,(0)) . (173) 
-OO 

Since the fields are independent and assumed to obey a free evolution, one finds: 

(<I>.(x)$,(y)) ^ TT[<^,{x)<^,{y)pp] = 5,, Gpix - y) , (174) 

where Gfj{x — y) is the standard li-dimensional Wightmann function for a single 
massless relativistic scalar field in a state at inverse temperature /3.^'' With the 
proper ie prescription, it can be written as: 

Gp{x)^ j ^-^0(fcO)<5(fc2)[(i+^^(^O))e-*-+AA;3(A;°)e*-]e-^'=\ (175) 
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where 



Its Fourier transform can be readily evaluated: 



/ + 00 
-OO 



(47r)Vr 



V 2 J 



C 1 



(176) 



(177) 



Remarks 4.1 



(1) As mentioned in Remark 3.7.3, the convergence in the weak coupling limit of the 
evolution equation for the reduced density matrix p{t) to the Markovian limit 
(172) stems from the asymptotic behavior of the two-point function Gfj{t) = 
Gjj{t,Q) for large t (see, (135)). More precisely, the combination \G i3{t)\{l + tY 
should be integrable on the positive half real line, for some r? > 0. In the case of 
a massless field, the integrals in (175) can be expressed in terms of the standard 
Poly-Gamma function, V^^H-^) = c^""^^ logT{z) / dz''+'^; explicitly, one finds:'^^ 

Since ^ 1/^" for \A ^ Arg(z) < tt, Gpit) behaves like l/t^~^ for 

large t in d > 4, while for d = 4 it falls off as due to a cancellation arising 
from taking the real part of V'*'^'' in (178). Therefore, the conditions for the 
existence of the weak coupling limit are assured, at infinity by the fall off of 
Gffit), provided > 4, and at zero by the ie prescription.' Notice that in two 
dimensions (178) gives a logarithmic behaviour for the two point function, so 
that in general the weak coupling limit is not guaranteed to exist. The limit is 
however well-defined for an atom in an optical fiber as discussed in Example 4.1. 
There, the field operators $ that couple with the atom variables involve the time 
derivative of a two-dimensional, free, massless scalar field; as a consequence, the 
relevant environment correlation ($(t)$(0)) actually behaves as d^Gf){t), and 
integrability at infinity is still assured. 

(2) The ie prescription involved in the definition of the Wightmann function (175) 
originates from causality requirements. In the present setting, it can be alterna- 
tively interpreted as an infinitesimal remnant of the size of the atom.^^ Indeed, 
assume that the atom has a spatial extension, with a profile given by a func- 
tion f{x). Since the atom-field interaction takes now place on the whole region 



G0{t) 



UTTi^rr^lK-/?) 



\d-2 



Tie 



In four-dimensional spacetime, the asymptotic behaviour for large space-hke separation of the 
two-point Wightmann function has been studied in Ref.[98] in the case of arbitrary interacting 
fields. In general, G/3(t) is found to fall off at infinity as so that the weak coupling limit is 

assured to exists even in this very general setting. 
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occupied by the atom, the field operator entering the Hamiltonian H' in (162) 
is now smeared out over the atom size, 



^iit) = J d^-^x f{x) ^i{t, x) , (179) 

instead of being just $»(*, 0), i.e. evaluated at the atom position. As a conse- 
quence, the correlation function (174) involves the Fourier transform f{k) = 
J d'^^^x c^'^'^ f{x) of the shape function. As an illustration, assume the fields to 
be in the vacuum, zero-temperature state; then, (compare with (175)); 

{Mmm = J ^ |M)|^ e-^l^^l* . (180) 

Since at the end we want to model a point-like atom, we choose a shape function 
approximating a (d — l)-dimensional (5- function: 



i=l 



[ix.r + {s/2y 



in this way, the atom is viewed as d — 1-dimensional box with size s. The 
Fourier transform of (181) can be easily computed: f{k ) = e~l'^l^/^; inserting 
this result in (180), one readily recovers the standard is prescription for the 
field correlations. □ 



As discussed in Section 3.3.2, the explicit expression of the Kossakowski matrix 
Cij involves a sum over the differences of the energy levels of the system Hamiltonian 
Hs, here just ±u)/2; it corresponds to the decomposition of the system operators 
CTj along the two energy eigenprojectors 

P± = i(l2±n-a) , (182) 

as indicated in (127), namely: 

CTi(O) = P+aiP+ + P-GiP- , a(±) = P±(JiP^ . (183) 

The 2x2 auxiliary matrices cri(O), (Ji{+) and <Ji{—) can be further decomposed 
along the Pauli matrices themselves, 

3 

= E ' for ^ = +, -, , (184) 

with the help of the three-dimensional tensor coeSicients 

V--"^ = m rij , ^ = ^ {^ij - 'f^i ^3 ± i^ijknk) ■ (185) 
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Recalling (173), one can then write (compare with (133)): 

3 

{=0,± kd=l 

= GpiO) + Gpiiu) ' + Gfji-u;) ^5+' , (187) 

where the second line follows from the property: X]fc=i '^'h ''i^kj ^ '^ij ■ ^'^^ ^^"^ 
easily check that, being the sum of three positive matrices, the matrix Cij is also 
positive, so that the dynamical semigroup generated by (172) consists of completely 
positive maps. As discussed before, this is to be expected, since (187) is the result 
of the weak coupling limit, in the mathematically well-defined formalism developed 
in Section 3.3.2. On the other hand, let us remark once more, that direct use of 
the standard second order perturbative approximation (as adopted for example in 
Ref.[100]) leads to physically inconsistent results, giving a finite time-evolution for 
p(t) that in general does not preserve the positivity of probabilities. 

Example 4.3 As an explicit illustration of this general framework, it is instructive 
to reconsider the simple system discussed in Example 3.3, and analyze its dynamics 
in the long-time, Markovian regime using the techniques just presented. The model 
describes the behaviour of a two-level atom immersed in a one-dimensional boson 
gas at temperature /3~^, that in the thermodynamical limit can be fully described 
in terms of a two-dimensional, free, massless scalar field. 

Positioning the atom at the origin, the total Hamiltonian in (81) can then be 
conveniently rewritten as: 
O 

H = -cT3®\4,^-ls®H^ + \(Ti®dt(l){Q) , (188) 

so that we are precisely in the condition of Eq. (161) and (162) with n = (0, 0, 1) and 
$1 = $2 = 0, <I>3 = dt(t>- The Fourier transform of the bath correlation functions, 
G/3(t) = Tr [9((^(t) dt(t>{^) P/s] , can be easily evaluated: 

through (186), one then obtains the expression for the Kossakowski matrix Cij. It 
turns out that only the entry C33 is actually nonvanishing. As a consequence, the 
Markovian master equation describing the reduced dynamics of the two-level system 
takes, in the interaction picture, the particularly simple form: 

dtPit) = ^[(73 Pit) as- Pit)] . (190) 

In the basis of eigenvectors of (73, CTsji) = (— l)*"*"^ |z), i = 0,1, one then recovers 
the exponential damping of the off-diagonal entries of pit) already discussed in 
Example 3.3: 

{i\pit)\z) = {M0)\i) , 

(l|p(i)|0)=e-2^^*/^l|p(0)|0). (191) 
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Although in a simplified settings, this result constitutes a direct evidence that the 
Markovian master equation (190), obtained using the mathematically well defined 
weak coupling limit procedure, correctly reproduces the long time behaviour of the 
exact reduced dynamics. □ 

The Lamb shift correction H^^^ to the effective Hamiltonian i?eff = Hs + 
can be similarly computed. It involves the Hilbcrt transform of the correlations 
(173), as indicated in (131) and (132). In practice, recalling (174), it can be expressed 
in terms of the following integral transform of the scalar Wightmann function 

^MC) = ^ r dz^. (192) 
ITT J_^ z-C, 



Explicitly, one finds: 



„(2) _ 1 ^ 



3 



Remark 4.2 This result deserves a closer look. Recalling (177), the definition of 
/C^(C) in (192) can be split as: 



^MC) = 



i(4.)^r(^) 



V / d^ 



d-3 



(194) 
(195) 



into a vacuum and a temperature-dependent piece. Although not expressible in 
terms of elementary functions, the temperature dependent second term is a finite, 
odd function of C„ vanishing as [3 becomes large, i.e. in the limit of zero temperature. 
The first contribution in (194) is however divergent for d > 3. As a consequence, 
despite some cancellations that occur in (193) (see below), the Lamb contribution 

(2) 

turns out in general to be infinite, and its definition requires the introduction 
of a suitable cutoff and a renormalization procedure. 

This is a well known fact, and has nothing to do with the weak coupling assump- 
tions used in deriving the master equation. Rather, the appearance of the diver- 
gences is due to the non-relativistic treatment of the two-level atom, while any sen- 
sible calculation of energy shifts would require quantum field theory techniques. 

In our quantum mechanical setting, the procedure needed to make the Lamb 
contribution well defined is therefore clear: perform a suitable temperature 
independent subtraction, so that the expression in (193) reproduces the correct 
quantum field theory result, obtained by considering the field in its vacuum state. □ 



60 Benatti, Floreanini 



4.2. Single Atom Dynamics and Decoherence 



From now on we shall work in ordinary four-dimensional Minkowski spacetime and 
therefore specialize d = 4. In this case, the Kossakowski matrix Cij in (187) takes 
the general form 

Cij = A 6ij - iB Cijk Uk + Cui Hj , (196) 

where the quantities A, B and C depend on the system frequency w and the inverse 
temperature /?: 

uj ri + e"^'^" 



A = 


1 

2 


B = 


1 
2 


C = 


1 



1 



UJ 
4-77 



UJ 

47r 



2 



1 + e- 



1 - e-P"^ 



On the other hand, the effective Hamiltonian iJgff can be rewritten as 

UJ 

Heff = - n-a , 
in terms of a redefined frequency: 

uj = u} + i [/C/3(— w) — /C/3(c<;)] . 



(197) 
(198) 
(199) 

(200) 
(201) 



Remark 4.3 As explained above, a suitable temperature independent subtrac- 
tion has been impUcitly included in the definition of the combination /C(— w) — /C(w), 
which otherwise would have been logarithmically divergent (compare with the ex- 
pression in (194)). It comes from a "mass" renormalization due to virtual vacuum 
effects of the full relativistic theory,^°^ and can not be accounted for by a mere 
frequency shift, as the expression (201) might erroneously suggest. □ 

In order to discuss explicitly the properties of the solutions of the master equa- 
tion (172), it is convenient to decompose the density matrix p in terms of the Pauli 
matrices, as introduced in Example 3.1: 

P=^(l2 + p-f?) . (202) 

Then, as shown there, the evolution equation (172) can be rewritten as a 
Schrodinger-like equation for the coherence (Bloch) vector \p{t)) of components 
(l,pi(i),P2(t),P3(t)): 

l\p{t)) = -2{n + V)\p{t)) . (203) 

The 4x4 matrices H and V correspond to the Hamiltonian and dissipative contribu- 
tions, respectively; they can be parametrized as in Example 3.1, Eqs.(53) and (54), 
in terms of the quantities Ui = {uj/2) n,, i = 1,2, 3, coming from i?eff) = — 4Bni, 
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V = —ABn2, w = —ABus, coming from the imaginary part of the Kossakowski 
matrix (196), and 

a = 2A + C{nl + nl) , b=-Cnin2, (204) 
a = 2A + C{nl + nl) , c=-Cnin3, (205) 
7 = 2A + C(nf +ni) , (3 = -Cn2n3, (206) 

coming from its real part. The 3x3, symmetric submatrix V^^^ defined in (56) 

with these last six parameters as entries, is manifestly positive; its eigenvalues are 
in fact Ai = 2A and X± =2A + C. As a consequence, for large times \p{t)) reaches 
a unique equilibrium state |/Oeg)-^°^ This asymptotic state is determined by the 
condition {H + 'D)\peq) = 0, and therefore is obtained by inverting V^^^: 

B 1 — e"''" 

\peq) = {l,-Rni,-Rn2,-Rn3] , R=a^iT^' ^^'^^^ 

Inserting these components in the expansion (202), one finds that the asymptotic 
density matrix is a thermal state in equilibrium with the field environment at 
the latter temperature T = 



s 



= T^^^ ■ (208) 

This asymptotic thermalization phenomenon is just one of the many effects that 
the dissipative dynamics generated by (203) produces on the two-level atom. More 
details on the thermalization process can be retrieved by studying the behaviour 
of the atom state p{t) for finite times; it is formally given by the exponentiation of 
(203): 

\p{t)) = e-2(«+^)* |p(0)) . (209) 

Recalling the definition (4), the evolution in time of any relevant atom observable, 
a Hermitian 2x2 matrix X, can thus be explicitly computed: 

{X{t))=Tv[Xp{t)] . (210) 

When the observable X represents itself an admissible atom state pf, the mean 
value (210) coincides with the probability Vi^f{t) that the evolved atom density 
matrix p{t), initially in p(0) = pi, be found in such a state at time t. Using (209), 
this probability can be computed explicitly: 

ViMt) = ^{l - (pf • n) (l - e-4^*^ - 



1 + c-l^"^ 



+ e-4^* (pi-n) (pen) 



+ e 2(2A+c)f ^ . _ . . cosuit - n • (pi X pf) sinCotj | . 

(211) 

In the above expression, we have expanded the density matrices pi, pf as in (202) 
and used the notations (pi • pf) and (pi x pf) (and similarly with n) to represent 
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scalar and vector products of their corresponding three-dimensional coherence vec- 
tors. As expected, this expression contains exponentially decaying factors involving 
the eigenvalues Ai and X± defined earlier, modulated by oscillating terms in the 
effective frequency Co. In other terms, the two-level atom immersed in scalar fields 
at temperature 1//3 is subjected to decoherence and dissipation, both effects being 
governed by the Planckian factors appearing in the Kossakowski matrix (196)-(199). 

When Pi = —n and pf = ft, the density matrices pi, pf represent the ground and 
excited states of the system Hamiltonian Hs in (161). In this case, the expression 
in (211) simplifies, 

^i^f W = (i - «"'^') ' (212) 

giving the probability for a spontaneous transition of the atom from the ground 
state to its excited state. Note that this phenomenon of spontaneous excitation 
disappears as /3 — > oo, i.e. when the environment is at zero temperature. 

Remark 4.4 Although the behaviour of 'Pi^{{t) in (211) and (212) is in principle 
experimentally observable through the use of suitable interferometric devices, in the 
study of the behaviour of atoms in optical cavities one often limits the discussion to 
the spontaneous excitation rate Fi^f, the probability per unit time of the transition 
i — > f, in the limit of an infinitely slow switching on and off of the atom-field 
interaction. In our formalism, its expression can be easily obtained by taking the 
time derivative of Vi^f{t) at t = 0; in the case of (212), one then finds 

Ti^f = - -^J—r . (213) 

□ 



4.3. Two Atoms: Environment Induced Entanglement Generation 

We have thus far analyzed the dynamics of an atom immersed in a bath of scalar 
masslcss fields at finite temperature, and found that, through decoherence effects, 
the atom state is driven towards a purely thermal equilibrium state, characterized 
by the same field temperature. When the subsystem in interaction with the field en- 
vironment consists of two, non-interacting two-level atoms one thus expects similar 
mixing-enhancing phenomena to occur, leading in particular to loss of the mutual 
quantum correlation (entanglement) that might have been present at the beginning. 

However, even though not directly coupled, the external field through which 
the two atoms move may provide an indirect interaction between them, and thus 
a means to entangle them. Indeed, entanglement generation through the action 
of an external heat bath has been shown to occur in certain circumstances; it is 
therefore of physical interest to investigate the same issue also in the present physical 
situation. 
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We shall therefore start by considering a system composed by two, equal two- 
level atoms, that start interacting with the external fields at time t = 0. Being 
independent, without direct mutual interaction, their internal dynamics can again 
be taken to be described by the generic Hamiltonian (161). Then, the total two- 
system Hamiltonian Hs is now the sum of two terms: 

3 3 

Hs = + Hf , H^p = I ^ „.(^. ^ 1) , iff = I ^ ^ . (214) 

i=l i=l 

Similarly, being immersed in the same set of free fields and within the weak coupling 
hypothesis, the atom- field interaction Hamiltonian can be most simply assumed to 
be the generalization of that in (162): 

3 

i?' = ^[(ai(g)l)(g)$i(a;) + (l(g)ai)(g)*i(a;) . (215) 

On the other hand, the Hamiltonian describing the environment dynamics remains 
that of a collection of free, independent scalar fields $8 and ^'j. 

Remark 4.5 In writing the interaction terms in (215), we have assumed the two 
atoms to interact with the field bath through two different field operators, $ and 
4', respectively; as wc shall sec, this allows discussing entanglement generation in 
a more general setting. For two equal atoms, the identification $ = \1/ is however 
physically justified, and we shall adopt it in most of the following analysis. On the 
other hand, since the two atoms are assumed to be pointlike, with infinitesimal size, 
there is no restriction in positioning them at the same point x, and the coupling in 
H' precisely reflects this simplified assumption.™ □ 

The derivation of the appropriate master equation describing the dynamics of 
the two atoms proceeds as in the case of a single atom, discussed in the previous 
Section. One starts from the Liouville-von Neumann equation (171) generating the 
time-evolution of the state /Otot(O) of the total system {atoms -|- external fields}, 
and then traces over the fields degrees of freedom, assuming a factorized initial 
state /9tot(0) = p(0) <8i P/3. In the weak coupling limit, the two-atom system state 
p{t), now a 4 X 4 density matrix, is seen evolving in time according to a quantum 
dynamical semigroup of completely positive maps, generated by a master equation 
in Kossakowski-Lindblad form: 

^ = -i[n,s,p{t)]+mt)]- (216) 

The unitary term depends on an effective Hamiltonian which is the sum of Hs 
in (214) and suitable Lamb contributions: HeS = H^s + H^s + The first 



""The case of spacially separated atoms can also be similarly treated, adopting the same open 
system techniques; for details, see Ref.[103]. 
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two represent single system contributions; they can be written as in (214), but in 
general with different terms for the two atoms: 

3 3 



The third piece is a field-generated direct two-atom coupling term, that in general 
can be written as: 



(218) 



As in the single atom case, a suitable temperature-independent renormalization 
procedure needs to be implemented in order to make all these contributions well 
defined. 

The dissipative contribution D[p(t)] can be written as in (49), 

6 



(219) 



with the help of the matrices = era (g) 1 for a = 1,2,3, = 1 (8) aa-3 for 
a = 4,5, 6. The Kossakowski matrix C^/j is now a positive 6x6 matrix, that can 
be represented as 



C = 



A B 
c 



(220) 



by means of the 3x3 matrices A = A\C = and B. In terms of this decomposition, 
can be rewritten in the following more explicit form: 

3 / 



»[p]=E 

i,j=i \ 

+Bij 



+b: 



(ctj l)p{cTi O 1) - ^|(a-i(Tj O 1) , p| 

(I0f7j)p(l0a-i)- ^|(l(g)C7iCrj), p| 
{(Tj O l)p(l®CTi) - ^{(cTj ®f7j), p| 

(1 ® (Tj)p((7j (g) 1) - ^|(o-i (g)aj) , p| 



(221) 



The structure of the various contributions reveals their direct physical meaning. 
Indeed, the first two contributions arc dissipative terms that affect the first, respec- 
tively the second, atom in absence of the other. On the contrary, the last two pieces 
represent the way in which the noise generated by the external fields may correlate 
the two, otherwise independent, atoms; this effect is present only if the matrix B is 
different from zero. 
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Remark 4.6 As in the single atom case, the environment contributions to Tieff 
and the entries of the matrix Cap in (220) can be expressed in terms of the Hilbert, 

respectively the Fourier transforms of the thermal Wightmann functions. In partic- 

f 121 

ular, the matrices ' in (218) and Bij in (221) do not vanish only if the bath state 
P/3 correlates bath-operators coupled to different atoms, that is, if the expectations 
Tr[/)^ ^i{t) ^j(O)] are nonzero. Only in this case, entanglement has a chance to be 
created by the action of the bath. Indeed, if -fff^ ^'^ = and Bij = 0, the two atoms 
evolve independently and initially separable states may become more mixed, but 
certainly not entangled. □ 

In order to investigate whether the thermal bath made of free fields can actually 
entangle two independent atoms, one can use the partial transposition criterion, 
that was introduced in Theorem 2.3. In fact, we are dealing here with a couple of 
two-level systems, and therefore Theorem 2.4 applies. In more precise terms, the 
environment is not able to create entanglement if and only if the operation of partial 
transposition preserves the positivity of the state p{t) for all times. 

From the operational point of view, one then prepares the two subsystems at 
t = in a separable state: 

P(0) = |V)(^|®|V)(^I ; (222) 

without loss of generality. p{0) can be taken to be pure: indeed, if the environment 
is unable to create entanglement out of pure states, it will certainly not entangle 
their mixtures. Then, one acts with the operation of partial transposition (over the 
second factor, for sake of definiteness) on both sides of (216). One thus obtains a 
master equation for the matrix p{t), the partially transposed p{t), that can be cast 
in the following form:^"^ 

dtp{t) = -i [Ues , Pit)] + B[p{t)] ; (223) 

here, Hcff is a new Hamiltonian to which both the unitary and the dissipative term 
in (216) contribute: 

3 3 

Weff = i^i ® 1) + E ^^'^ (1 ® 

i=\ ij=l 

3 

+ Y,^m{B-E)^^ia,^aj) , (224) 

ij=l 

where E is the diagonal 3x3 matrix given hy: E = diag(— 1, 1,-1). The additional 
piece D[ • ] is of the form (219), but with a new matrix C ^ £ ■ C ■ £, where 

C-( ^ ne{B) + ^H(^-)\ 



I3 
E 



(226) 
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and ^ denotes full matrix transposition, while iJ^^^) is the coefficient matrix in (218). 

Although p{t) evolves according to a master equation formally of Kossakowski- 
Lindblad form, the new coefficient matrix C need not be positive. As a consequence, 

the timc-cvolution generated by (223) may result to be neither completely positive, 
nor positive and need not preserve the positivity of the initial state p{0) = p{0). 

Remark 4.7 Notice that both the Hamiltonian and the dissipative terms of the 
original master equation (216) contribute to the piece ©[-J in (223), the only one 
that can possibly produce negative eigenvalues. In particular, this makes more trans- 
parent the physical mechanism according to which a direct Hamiltonian coupling 

('12) 

H^^ among the two systems can induce entanglement (as studied for instance in 

('12) 

Refs.[104-109]): on p{t), H^^' "acts" as a generic dissipative contribution, which in 
general does not preserve positivity. □ 

In order to check the presence of negative eigenvalues in p{t), instead of exam- 
ining the full equation (223) it is convenient to study the quantity 

Q{t) = {x\~p{t)\x) , (227) 

where x is any 4-dimensional vector. Assume that an initial separable state p has 

indeed developed a negative eigenvalue crossing the zero value at time t* . Then, 
there exists a vector state |x) such that Q{t*) = 0, Q{t) > Qioit <f* and Q{t) < 
for t > t*. The sign of entanglement creation may thus be given by a negative first 
derivative of Q{t) at t = t* . Moreover, by assumption, the state p{t*) is separable. 
Without loss of generality, one can set i* = and, as already remarked, restrict the 
attention to factorized pure initial states. 

In other words, the two atoms, initially prepared in a state p(0) = p(0) as in 
(222), will become entangled by the noisy dynamics induced by their independent 
interaction with the bath if there exists a suitable vector |x), such that: 

1) Q(0) = and 

2) dtQ{0) < 0. 

Note that the vector |x) needs to be chosen entangled, since otherwise Q{t) is never 
negative. 

Remark 4.8 The criterion just stated is clearly sufficient for entanglement creation. 
In fact, when 9(Q(0) > for all choices of the initial state p(0) and probe vector |x), 
entanglement can not be generated by the environment, since p remains positive. 
However, the treatment of the case dtQ{Q) = requires special care: in order to 
check entanglement creation, higher order derivatives of Q, possibly with a time 
dependent |x), need to be examined. □ 

A more manageable test of entanglement creation, valid for any probe vector 
Ix), can be obtained by a suitable manipulation of the expression 9tQ(0). In the 
two-dimensional Hilbert spaces pertaining to the two atoms, consider first the or- 
thonormal basis {|<^), |<^)}, {l^i"), |'*/^)}, obtained by augmenting with the two states 
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\(p) and 1-0) the ones that define p{{)) in (222). They can be both unitarily rotated 
to the standard basis {|— ), |+)} of eigenvectors of CT3: 

\^) = U\-) 1^) = U\+) , 

\i^)=V\-) m=V\+). (228) 

Similarly, the unitary transformations U and V induce orthogonal transformations 
U and V, respectively, on the Pauli matrices: 

3 3 

U^aiU = Y^U^jaj , V^'aiV = ^ VijCFj . (229) 

Direct computation then shows that 9(Q(0) can be written as a quadratic form in 
the components of the probe vector |x). As a consequence, vectors |x) exist making 
this form negative, i.e. 9fQ(0) < 0, if and only if its corresponding discriminant is 
negative; explicitly: 

{u\A\u){v\C'^\v) < \ {u\{ne{B) + iH^^'^'^)\v)\^ . (230) 

The three-dimensional vectors \u) and \ v) contain the information about the starting 
factorized state (222): their components can be in fact expressed as: 

3 3 

«, = ^Z^,^.(+|a,|-) , Vi = Y,Vij{-\cJ,\+) . (231) 

Therefore, the external quantum fields will be able to entangle the two atoms evolv- 
ing with the Markovian dynamics generated by (216) and characterized by the Kos- 
sakowski matrix (220), if there exists an initial state \'^){'^\ \'tp){ip\, or equivalently 
orthogonal transformations U and V, for which the inequality (230) is satisfied. 

The test in (230) is very general and can be applied to all situations in which two 
independent subsystems are immersed in a common bath. It can be satisfied only 
if the coupling between the two subsystems induced by the environment through 
the mixed correlations Tr[p^$i^'j], encoded in the coefficients of B and H^^'^\ are 
sufficiently strong with respect to the remaining contributions. 

As already remarked, in the specific case of two atoms in interaction with the 
same set of external fields, it is physically reasonable to take = so that 
the just mentioned mixed, environment generated couplings result maximal and 
equal to the diagonal ones, i.e. Aij = Bij = Cij.^ Furthermore, in such a case, 
recalling (174), all field correlations become proportional to the standard, scalar 
field temperature Wightmann function Gf}{x) introduced in (175); as a consequence, 
the three matrices A, B and C in (220) become all equal to the Kossakowski matrix 
(196) pertaining to a single atom: 

Aij = Bij = Cij = A Sij ~ iB eijk Uk + C rn Uj . (232) 

"This particulcir choice for the Kossakowski matrix (220) is of relevance also in phenomeno- 
logical applications; for instance, it is adopted in the analysis of the phenomenon of resonance 
fluorescence.^^"'^* 
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In such a situation, also the Hamiltonian pieces (217) and (218) simplify: H^J> and 

(2) 

H^g can be written exactly as the original system Hamiltonian (214), with the 
frequency u) replaced by the renormalized one u) given in (201), while the direct 
two-atom coupling term takes the form 

^iff ^ = E {[/C/3(a;) +/CM-a;)] Sij 

+ [2/C/3(0) - ICf3{aj) - ICf3{-Lo)] niiij^ at (g) CTj , (233) 

with /C^(C) as in (192). 

Remark 4.9 Recall that this function can be split as in (194) into a vacuum and 
temperature dependent piece. Since, as observed there, the temperature dependent 

f 12) 

contribution to K. (}{(,) is odd in C,, one deduces that H^^ ' does not actually involves 
T = 1//3: it is the same Lamb term that would have been generated in the case of a 
two-atom system in the vacuum. Being interested in temperature induced phenom- 
ena, we shall not consider this vacuum generated term any further and concentrate 
the attention on the effects produced by the dissipative contribution D[p] in (219). □ 



Let us then consider again the test for entanglement creation given in (230). 
Because of condition (232), it now involves just the hermitian 3x3 matrix Aij. 
Its expression can be further simplified by choosing Ui = Vi in (231); recalling the 
definitions (228), (229), this in turns implies = \(f)) (in other terms, if in the 
initial state p{0) we choose |^) = |— ), then l'^) must be taken to be |+)). In this 
case, (230) simply becomes: 

\{u\Im{A)\u)\'^ > . (234) 

As long as Aij is not real, i.e. the parameter B in (232) is nonvanishing, this 
condition is satisfied for every |m) outside the null eigenspace of Jm(^). Take for 
instance the initial state p(0) = |— ) (— (+|, so that the three-dimensional vector 
\u) has components Ui = {1, — i,0}; from (232) one easily finds: | (w|Ito(^)|u) | = 
{Bn^Y' , which is in general non vanishing. 

We can thus conclude that entanglement between the two atoms is indeed gener- 
ated through the weak coupling with the external quantum fields. This happens at 
the beginning of the time-evolution, as soon as t > 0. Note however that the test in 
(234) is unable to determine the fate of this quantum correlations, as time becomes 
large. In order to discuss asymptotic entanglement, one has to analyze directly the 
structure of the dynamics generated by the master equation in (216). 
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4.4. Two Atom Reduced Dynamics: Entanglement Enhancement 

As in the case of the single atom evolution, it is convenient to decompose the 4x4 
density matrix p{t) describing the state of the two atoms along the Pauli matrices: 



^(*) = 4 



3 3 3 

1 (g) 1 + ^ PQi{t) 1 (g) (Tj + '^pio{t) (Tj 1 + ^ Pij{t) (Ti (Jj 
i—1 i—1 ij — 1 



, (235) 



where the components poi{t), Pio{t), Pij{t) are all real. Substitution of this expansion 
in the master equation (216) allows deriving the corresponding evolution equations 
for the above components of p{t). As explained in Remark 4.9, wc shall ignore the 
Hamiltonian piece since it can not give rise to temperature dependent entanglement 
phenomena, and concentrate on the study of the effects induced by the dissipative 
part in (221), with the elements of the Kossakowski matrix as in (232). 

Let us first observe that when the three submatrices A, B and C in (220) are 
all equal, the the form of the dissipative contribution in (221) simplifies so that the 
evolution equation can be rewritten as 



dpit) 
dt 



3 r ^ 

= ]D>[p{t)] = ^ Aij E,- p{t) - -{SiEj , p{t)} 



in terms of the following symmetrized two-system operators 

= O-j (g) 1 + 1 (g) (Tj 2 = 1,2,3. 



(236) 



(237) 



One easily checks that these operators obey the same su{2) Lie algebra of the Pauli 
matrices; further, together with 



(Ti C7i + C7i (7,; 



ij = 1,2,3 , 



(238) 



they form a closed algebra under matrix multiplication. For later reference, we 
give below its explicit expression: it can be easily obtained using aiaj = Sij + 

* Z^fe=l ^ijk Cfe- 

3 

Ei Ej = 2 1 1 + z ^ Eijk Efe + Sij , 
fe=i 

3 3 

Sik Ej + Sjk Ej 

+ i ^ Siki Sij + i ^ Ejki Sii , 
1=1 1=1 

3 3 

Sik Ej + Sjk ^i — i ^ Siki Sij — i ^ Ejki Sii , 
1=1 1=1 

3 

Sij Ski = 2(i5ifc Sjl + 5ii Sjk)l 1 + « ^ (^SikSjlr + SjkSilr + SilEjkr + SjlSikr^ E^ 



Sij Efc 



Efc Sij — I 



- ^2 6ij 6ki — Sik Sjl — Sii Sjk^ S + 2^Sij Ski + Ski Sij^ 
-Sik Sjl — Sii Sjk — Sjk Sii — Sjl Sik I 



(239) 
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where S = X]r=i '^a- 

Inserting the exphcit expression for Aij given in (232) into (236) and using the 
decomposition (235), a straightforward but lengthy calculation allows to derive the 
following evolution equations for the components of p{t): 

3 

-2^ { [{2A + C) Sik - Cn.nk] Pok{t) - Buk p^k{t)\ ~ 2B{2 + r) n, , 



dpoijt) 
dt 

dpio{t) 



fe=i 

3 



dt 



-2J2{ [(2^ + C) Sik - Criirik] Pko{t) - Buk Pk^{t)} - 2B{2 + r) n, 
fe=i 

= -4[{2A + C) pij{t) + {A + C) pjiit) - ((A + C) 5ij - CmUj), 
-AB[ni poj{t) + rij pio{t)] - 2B[ni pjo{t) + rij poiit)] 

3 

+'^^ [b Sij nk[pko{t) + pok{t)] +Cnink[pkj{t) +2pjk{t)] 



fe=i 

3 



+Cnjnk[pikit) +2pki{t)]^ - 4C ^ UkUi pki{t) . (240) 



feJ=i 



In these formulae, the quantity r = Y^^=i Pa represents the trace of Pij; it is a con- 
stant of motion, as easily seen by taking the trace of both sides of the last equation 
above. Despite this, the value of t can not be chosen arbitrarily; the requirement of 
positivity of the initial density matrix p{0) readily implies: — 3 < r < 1. 

The system of first order differential equations in (240) naturally splits into two 
independent sets, involving the symmetric, p(oi) = Poi + Pio, P{ij) = Pij-^Pjii a-nd an- 
tisymmetric, p[oi] = poi — Pio, P[ij] = Pij — Pji, variables. By examining the structure 
of the two sets of differential equations, one can conclude that the antisymmetric 
variables involve exponentially decaying factors, so that they vanish for large times. 
Then, using the definitions (237) and (238), the study of the equilibrium states p 
of the evolution equation (236) can be limited to density matrices of the form 



^ 1 
^=4 



Pij Sij 



(241) 



with Pij = Pji. 

The approach to equilibrium of semigroups whose generator is of the generic 
Kossakowski-Lindblad form has been studied in general and some rigorous math- 
ematical results are available. ^^^'^ We shall present such results by adapting them 
to the case of the evolution generated by the equation (236). 

First of all, one notices that in the case of a finite dimensional Hilbcrt space, 
there always exists at least one stationary state po: this can be understood by 
recalling that in finite dimensions the ergodic average of the action of a completely 
positive one-parameter semigroup on any initial state always exists; the result is 
clearly a stationary state. 
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Let us now introduce the operators Vi = ^1/^ (recall that A is non- 

negative), so that the r.h.s. of (236) can be rewritten in diagonal form: 

3 p ^ 

»[P] = E ^0 - ^{^1^0 ^ P) ■ (242) 

i,j = l 

When the set M formed by all operators that commute with the linear span of 
{Vi, Vj, i = 1,2,3} contains only the identity, one can show that the stationary 
state po is unique, and of maximal rank. On the other hand, when there are several 
stationary states, they can be generated in a canonical way from a po with maximal 
rank using the elements of the set A4 . 

In the case of the evolution equation (236), A4 contains the operator S = 
Y^^=i^ii' besides the identity; indeed, with the help of the algebraic relations in 
(239), one immediately finds: [S, Sj] = 0. Out of these two elements of M, one can 
now construct two mutually orthogonal projection operators:" 



p = i 

4 



1 - 

2 



g = 1 - P . (243) 



Then, one can show that any given initial state p(0) will be mapped by the evolution 
(236) into the following equilibrium state: 

That this state is indeed stationary can be easily proven by recalling that P and Q 
commute with Sj, i= 1, 2, 3, and thus with the Vj as well; therefore, P[/3] = 0, for 

any p{0), as a consequence of ID>[po] = 0. 

The problem of finding all invariant states of the dynamics (236) is then reduced 
to that of identifying a stationary state po with all eigenvalues nonzero. Although 
in principle this amounts to solving a linear algebraic equation, in practice it can 
be rather difficult for general master equations of the form (172). Nevertheless, in 
the case at hand, the problem can be explicitly solved, yielding: 

Po = ^{U - Rn ■ a'^ ^ i(l2-Pn-a) , (245) 

where R = B/A is the temperature dependent ratio already introduced in (207): 
note that < i? < 1, where the two boundary values correspond to the infinite and 
zero temperature limits, respectively. 

Inserting this result in the expression (244) allows deriving the expression of 
the set of all equilibrium states of the dynamics (236); as expected, they take the 
symmetric form of (241), with the nonvanishing components given by: 



°One easily checks that P is the projection operator Q_ on one of the maximally entangled Bell 
states introduced in Example 2.3. 
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1 



(r-iJ^) 5ij+R^{T + 3) mrij 



(246) 



Pij = 



2(3 + i?2) 



These stationary density matrices depend on the initial condition p(0) only through 
the value of the parameter r, that as already mentioned is a constant of motion for 
the dynamics in (236). 

As we have explicitly shown before, an environment made of the quantum free 
fields is in general able to initially entangle two independent subsystems weakly 
interacting with it. Nevertheless, by examining the dynamics of a subsystem made 
of a single atom, we have also seen that such an environment produce dissipation 
and noise, leading to effects that generically counteract entanglement production. It 
is therefore remarkable to find the asymptotic state p in (246) to be still entangled. 

To explicitly show this, one can as before act with the operation of partial 
transposition on p to see whether negative eigenvalues are present. Alternatively, 
one can resort to concurrence as discussed in Section 2; indeed, the expressions 
in (246) are sufficiently simple to allow a direct evaluation. The use of concurrence 
has another advantage; it not only signals the presence of entanglement, it also 
provides a measure of it: its value ranges from zero, for separable states, to one, for 
fully entangled states, like the Bell states. 

As explained in Theorem 2.5, in order to determine the concurrence of any 4x4 
density matrix p representing the state of two atoms, one computes the eigenvalues 
of the auxiliary matrix p {02 ® 02) p* {02 ® 02), which turn out to be non- negative; 
their square roots A^u, /i = 1, 2, 3, 4, can be ordered decreasingly in value: Ai > A2 > 
A3 > A4. The concurrence of p is then defined to be: C[p] = max{Ai — A2 — A3 — A4, 0}. 

In the case of the asymptotic state p in (246), the above mentioned procedure 
gives: 



This expression is indeed nonvanishing, provided we start with an initial state p(0) 
for which 



The concurrence is therefore a linearly decreasing function of t, starting from its 
maximum C = 1 for r = — 3 and reaching zero at r = (Si?^ — 3)/(3 — R^)\ notice 
that this ratio is an admissible value for r, since it is always within the interval 
[—1,1] for the allowed values of R. 

This result is remarkable, since it implies that the dynamics in (240) not only can 
initially generate entanglement: it can continue to enhance it even in the asymptotic 
long time regime. 




(247) 
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Example 4.4 As initial state, consider the direct product of two pure states: 

P{0) = Px <Si py , ps=^(l2 + x- aj , pg=^(l2 + y-a^ , (249) 

where x and y are two unit vectors. In this case, one easily finds that t = x ■ y, so 
that, recalling (247) the asymptotic entanglement is maximized when x and y are 
coUinear and pointing in opposite directions. In this case on explicitly finds: 

C[p] = ^ , (250) 

provided R ^ 0. Notice that C[p] reaches its maximum value of 1/2 when R = 1, 

i.e. at zero temperature, while it vanishes when the temperature becomes infinitely 
large, i.e. R= 0. This has to be expected, since in this case the decoherence effects 
of the bath become dominant. □ 



Furthermore, one can easily check that the phenomenon of entanglement pro- 
duction takes place also when the initial state p(0) already has a non-vanishing 
concurrence. Let us consider the following initial state, 

p(O) = |l0l + (l-e)P, (251) 

which interpolates between the completely mixed separable state (see Example 2.3) 
and the projection P in (243). Provided e < 2/3, this state is entangled, with 
= 1 ~ 3£/2. As this initial state evolves into its corresponding asymptote p, 
the difference in concurrence turns out to be 

C[p]-C[p(0)] = ^^, (252) 

which is indeed non vanishing. By direct inspection, one sees that the state P above 
is a fixed point of the dynamics generated by (240) and therefore it coincides with 
its corresponding asymptotic state p. This is in agreement with the already observed 
fact that maximal concurrence require r = — 3. 

Remark 4.10 The formalism used to analyze the behaviour of stationary atoms 
in a thermal field bath can be adapted to study the dynamics of a uniformly accel- 
erating point-like detector through a vacuum scalar field. Indeed, the detector can 
be modeled as a two-level system, in weak interaction with the external, relativistic 
quantum field, which plays the role of environment. 

In the comoving frame, the detector is seen evolving with a master equation of 
the form (172) as if it were immersed in a thermal bath, with temperature pro- 
portional to its proper acceleration. This phenomenon is known in the literature 
as the Unruh effect. The discussion of the previous sections suggests that for an 
accelerating subsystem composed of two, non-interacting two-level atoms, quantum 
correlations should be generated between the two subsystems as a result of the ef- 
fective open system dynamics. A detailed treatment, similar to the one outlined in 
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the last Section, precisely confirms that in general the asymptotic density matrix 
describing the equilibrium state of the two atoms turns out to be entangled. □ 

In summary, we have seen that the entanglement generating properties of a bath 
of quantum fields at finite temperature can be ascertained through the study, with 
different techniques, of two separate time regimes in the evolution of two subsystems 

immersed in it. Entanglement production in the short time region is regulated by the 
properties of the generator of their subdynamics, while asymptotic entanglement is 
measured by the concurrence of the final equilibrium state. 

The same techniques can also be used to study entanglement production in more 
general environments, for which the two subsystems time-evolution is still generated 
by master equations of type (236), but with a generic matrix Aij, not necessarily of 
the form (232).^^^ These extended dynamics are useful in phenomenological appli- 
cations, in particular in quantum optics:^^'''^ the fact that also in these generalized 
situations the asymptotic entanglement results non-vanishing may give important 
feedback for the actual realization of elementary quantum computational devices. 
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